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r ~~ 1 Abstract 

Q_i This document contains a review on the quark model, prepared for lectures at the 

o | Niccolo Cabeo School at Ferrara in May 2012. It includes some historical aspects, 

the spectral properties of the 2-body and 3-body Schrodinger operators applied to 
mesons and baryons, the link between meson and baryon spectra, the role of flavour 
independence, and the speculations about stable or metastable multiquarks. The 
analogies between few-charge systems and few-quark bound states will be under- 
\Q lined. 
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1 Prelude: few charge systems in atomic physics 

Si parva licet componere magnis^ 

Virgil 

1.1 Introduction 

The spectrum of few-charge systems was among the very first applications of quantum 
physics. The Bohr-Sommerfeld rules explain the energy levels of the hydrogen atom 
(j>, e~), and can be easily extended to any (mf, m^) pair with arbitrary masses. The so- 
lution of the three-body problem for (mf , , ) turned out less obvious, and required 
the Schrodinger equation, and the associated variational methods. It also revealed some 
surprises, with the stability imposing constraints on the masses 

The stability of the positronium molecule (e + , e + , e~ , e~), or any similar system with 
equal masses (in the limit where annihilation is neglected) was suggested by Wheeler 
in 1945 (IJ and proved in 1947 |2j, but the (indirect) experimental evidence was pub- 
lished only in 2007 [3 |. This indicates how patient one should be when predicting novel 
structures. 

The quantum mechanics of few-charge systems was a great source of inspiration 
for building the quark model, in its minimal non-relativistic version. Amazingly, some 
techniques developed to extrapolate into flavour space the few-quark spectra with fla- 
vour-independent forces turned out useful to understand the stability patterns in atomic 
physics. I refer to El for a review about few-charge systems. Here, I will just stress a 
few points that are connected to the quark model. 



1.2 The atomic two-body problem 
1.2.1 Central potential 

Out of the Hamiltonian 

2^ + 2^--' »> 
2 rrii 2 mi r 12 

the centre of mass can be removed and one is left with the intrinsic Hamiltonian de- 
scribing the relative motion, which reads 

H=f-~, (2) 
2 p r 

where p is the reduced mass, given by p^ 1 = m{ 1 + raj 1 , and r = r 2 — r 1 and p = 
(p 2 — Pi)/2 are conjugate variables for the relative motion. 

The Hamiltonian H has seemingly two parameters, the reduced mass p and the 
strength e 2 , but it can be rescaled. It is sufficient to solve once for ever the universal 
spectral problem 

h = -A - r- 1 , (3) 



Hf it is allowed to compare small things with great 
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and to apply a simple factor 2 /i e 4 to the eigenvalues and (2 ji e 2 ) -1 to the distances to 
recover the actual energies and wave-functions of 

There is a similar universality for all harmonic oscillators p 2 /m + Kr 2 , and more 
generally for all power-law interactions |5j 

H=^- + ge(a)r\ (4) 

where g > 0, e is the sign function, to ensure that the force is attractive, and a is not too 
large if negative, otherwise the Hamiltonian would not be bounded from below. 

A common feature of the Coulomb and Hookes potentials in three dimensions, is 
that they support an infinite number of bound states, however weak is the strength of 
the potential. This contrasts with short-range potentials such as the Yukawa interaction 
— g exp(— br)/r of nuclear physics, which requires a minimal strength g to achieve bind- 
ing. Similarly, the effective potential between neutral atoms or the effective potential 
between two hadrons are of short-range nature, and thus do not necessarily support 
bound states even when they are attractive. 

The Coulomb Hamiltonian Q has a well-known spectrum, with a discrete part e n> i = 
— 1/[4 [n + £) 2 }, where n = 1, 2, ... is the radial number, and I = 0, 1 . . . the orbital 
momentum (n+£ is the principal quantum number). There is also a continuum for e > 0. 



1.2.2 Spin-dependent corrections 

The Coulomb interaction can be derived from Quantum ElectroDynamics (QED) in the 
non-relativistic limit, and in the approximation of the lowest order in the coupling con- 
stant, which corresponds to one-photon exchange. 

There are several interesting corrections, which have been probed successfully. In 
particular, the vector nature of the photon gives very characteristic spin corrections. 

For £ > 1, there are spin-orbit and tensor terms which contribute to the fine structure: 
levels with the same orbital momentum £ but different coupling of £ to the spins have 
slightly different energies. 

For £ = 0, there is spin-spin term, which splits the spin-singlet form the spin triplet 
states. It reads 

V ss = — ^V)^, (5) 

in its most simplified form, to be treated at first order in perturbation theory, using the 
wave functions generated by the central potential —e 2 /r. The short-range character, and 
the dependence on the masses should be kept in mind when an analogue will be pro- 
posed for the quark-antiquark interaction. This interaction is responsible for the shift 
between the ortho- and para-hydrogen states and the transition between them produces 
the famous 21 cm line, whose gravitational red-shift and Doppler shift gives valuable in- 
formation in astrophysics. The analogue will be measured at CERN for antihydrogen, 
to probe the matter-antimatter symmetry. 
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Note that in the case of the positronium atom, the hyperfine splitting also receives 
some contribution from the annihilation diagrams. 



1.3 Three-unit-charge ions 

The best known 3-body system is atomic physics is the neutral Helium atom, (a, e~, e~), 
but its stability is obvious since once the first electron is attached to the a particle, there 
remains some long-range attraction to bind the second electron. 

The binding of H = (p, e~,e~) is less obvious, as the second electron feels a neutral 
system at large distances. Unlike the case of the Helium atom, neither an effective- 
charge ansatz ip(r 2 ,r 3 ) = exp(— a(r 2 + r 3 )), nor a more general Hartree wave function 
f(r 2 ) f(r 3 ) suffices to demonstrate the binding variationally. But a more elaborate wave 
function does. For a review on two-electron systems, see, e.g., J6j. 

More interesting, perhaps, is the dependence upon the masses |4j: the molecular 
hydrogen ion H 2 + = (p,p, e~), is very stable, with a variety of excitations, the hydrogen 
ion H = (p, e~, e~) and the positronium ion Ps _ = (e + , e", e") are bound by a small 
margin, but the "protonium ion" (p,p, e^) is not stable. (Here p and p simply denotes 
the mass and the Coulomb charge; any hadronic interaction is neglected.) 

There is an increase of stability near the symmetry line m 2 = m 3 of (mf ,rri% ,mf). 
Schematically, the system has there two threshold configurations (m|,m 2 ) + m 3 , and 
(mf, mj) + m 2 , which are degenerate, and thus interfere optimally to achieve the best 
binding. 

This mass dependence is sometimes acute. For instance, a fictitious H with one 
electron heavier by about 10% would not be be stable. This is the case, e.g., for the 
(very) exotic (p, pT , tt~) ion. 



1.4 Three-body exotic ions 

It is sometimes believed that H has no excited bound state. This is both true and false. 
This is true if you define a bound state as lying below the lowest threshold, which is 
made of H = (p, e~) in the IS state and an isolated electron. However, there is a level 
with total orbital momentum L = 1 and parity +1, that cannot break into H(1S) + e~ 
without involving spin forces and radiative effects. For this state, the lowest threshold 
for spontaneous dissociation consists of H = (p, e ) in the 2P state and an electron. A 
bound state lying below this threshold exists for H = (p, e~, e ), but it exhibits even 
a more striking mass dependence when the constituents are modified \7\. Allowing a 
very small mass difference between the two electrons (less than 1%) or replacing the 
proton by a positron spoils the binding. On the other hand, if the masses are inverted, 
i.e., for the unnatural-parity state of H 2 + = (e~,p, p), a comfortable binding is observed. 
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1.5 Molecules with four unit charges 

The best known case corresponds to the hydrogen molecule H 2 = (p,p, e~,e~), bound 
well below the threshold for dissociation into two hydrogen atoms, with a variety of 
excitations. There are also variants with one or two protons replaced by an isotope. The 
most suited framework for studying H 2 is the Born-Oppenheimer approximation. For 
a given position of the protons, the electronic energy is computed, and interpreted as an 
effective proton-proton interaction supplementing their direct electrostatic repulsion. 
The Schrodinger equation is then solved for the two-body problem of the protons and 
generates the ground state and its radial and orbital excitations. Another set of levels 
correspond to the electron cloud being in an excited state. (Later the analogue will be 
the hybrid hadrons, with the gluon field linking the quarks being excited). Note that 
one hardly treats the hydrogen molecule as a "di-proton" bound to a "di-electron". 

In atomic physics and ab-initio chemistry, one usually starts from the limit where the 
proton and other nuclei are infinitely heavy, and then their finite mass can be treated as a 
correction, through the Hughes-Eckart terms or similar. In 1945, Wheeler JTJ addressed 
the somewhat opposite issue of a proton with the same mass as the electron, i.e., the 
problem of the stability of the positronium molecule, (e + , e + , e~ , e~) (in the limit where 
annihilation is neglected). In 1946, Ore performed an elaborate 4-body calculation, and 
concluded that such a configuration is hardly stable [8|. But the following year, Hyller- 
aas and the very same Ore published an elegant analytic proof of the stability of the 
positronium molecule pj. It took 60 years to obtain at last an indirect experimental 
evidence of the existence of this molecule (3). 

Meanwhile the dependence upon the masses was studied in some details. In Refs. |9 



12 [, it was shown how (M + , M + , m~ , mr ) gets more binding, as compared to its thresh- 
old, when the mass ratio departs from M/m = 1. On the other hand, it was observed 
that when allowing for two different masses in another way, (M + , m + , M~,m~), stabil- 
ity is lost for M/m > 2.2 (or, M/m < 1/2.2), because the molecule cannot compete any 



more with the lowest threshold (M+, M~) + (m+, m~) |13||14| . 

Many less symmetric configurations are bound, such as the positronium hydride 
PsH = (p, e + ,e~, e~), or any configurations in which two particles have the same charge 



and the same mass, that is to say, (a + , b + ,c , c ), whatever are the masses of a and b | 15 [. 
The lessons, to be kept in mind when arriving at the section on multiquarks, are: 

• the role of the masses is important. Look at (a + , b + ,c~ ,d~) vs. its lowest threshold 
supposed to be (a + , c~) + (b + ,d~). For both the 4-body system and its threshold, 
the potential energy r^ 1 has a cumulated strength g {j = — 2. So, there is 
no obvious reason why a molecule should lie lower that two atoms. Hence any 
favourable breaking of symmetry in the kinetic energy is welcome. Conversely, 
any unfavourable breaking can spoil the binding. 

• the four-body problem is delicate, even when the interaction is known perfectly, 

• one should be patient to see the experimental confirmation of theoretical predic- 
tions of exotic states. 
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2 A brief historical survey 



History is Philosophy teaching by examples^ 

Thucydid 

2.1 Prehistory 

There are several books relating the birth of particle physics, with entertaining anec- 
dotes. Segre, for instance, was an acute observer |16) . A very comprehensive review 
is given by Pais |17[. For the nuclear forces, see |18[. Books and reviews devoted to 



the quark model include 1 19-25). There are also reviews covering baryons fl26j-|28) or 



mesons |29 -31 1, in particular the synthesis by the "quarkonium working group"" |32p3[ . 

After the Rutherford experiment, indicating how compact is the atomic nucleus, and 
the discovery of the neutron by Chadwick, it became necessary to understand how the 
nucleus is built out of protons and neutrons. The mechanism of Yukawa: the exchange 
of a massive boson, turned out the be successful, with the discovery of the pion at Bristol 
in 1947. 

Among the theoretical activity stimulated by the Yukawa model, two points at least 
deserve attention. 



1. The mass of the Yukawa boson is constrained by the ratio of the 3-body to the 
2-body binding energy, as stressed by Thomas in a celebrated paper |34"| , who 
anticipated what is known today as "Borromean binding", and, more generally, 
"Efimov physics". 

2. At first sight, three potentials are needed to build the nucleus: proton-proton, 
proton-neutron and neutron-neutron, and it is natural to seek for some simplifi- 
cation. A tempting scenario is that where solely the proton-neutron interaction 
exists, but this is clearly contradicted by the data on proton scattering. What even- 
tually prevailed is isospin symmetry. With the proton and the neutron forming a 
doublet, there are only two potentials in the limit where isospin is conserved, one 
for total isospin 1 = 0, and another one for 1 = 1. 



2.2 Early hadrons 

The pion was discovered in 1947 and seen in three charge states, n + , n°, and ir~ which 
form an isospin triplet. Hence in 1947, we had 5 hadrons: 2 nucleons and 3 pions. 
But already two more were expected, as the existence of antimatter was (not so easily) 
predicted as a consequence of the Dirac equation, and the positron was discovered in 
cosmic rays by Anderson. The antiproton and the antineutron were anticipated as well. 
A dedicated accelerator was built at Berkeley, the Bevatron, and the antiproton was, 
indeed, discovered in 1955, and the antineutron shortly after. 



2 According to Michel Casevitz, the sentence is not by Thucydid, but a British commentator 
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With 7 hadrons, 2 nucleons, 2 antinucleons, and 3 pions, the world of hadronic 
physics would be reasonably sized, and one could envisage to work on the interaction 
among these few hadrons. However, several complications occurred almost simultane- 
ously 

First, the Yukawa picture of nuclear forces, though very efficient for the long-range 
part, faced difficulties at shorter distances. More attraction was needed, and also some 
spin-orbit component, that neither pion-exchange or iterated pion-exchange were able 
to provide. The work by Breit, among others, was remarkable |35[. For further details, 
see, e.g., the review |18|. An explicit scalar exchange (call it a or e) and an explicit 
vector exchange were needed. While the former was about isospin independent, and 
thus provided by the exchange of an isoscalar scalar meson, the latter was sought to be 
different in np vs. pp scattering, and thus called for both an isoscalar and an isovector 
vector meson. The to and p were thus predicted! 

Second, the interaction of pions with nucleons was shown to produce new parti- 
cles, nucleon resonances, in particular the A(1232), which has isospin 3/2, i.e., exists in 
four possible electric charges. Similarly, proton-nucleon or proton-nucleus scattering, 
or proton-antiproton annihilation were able to produce several new mesons, the ones 
desired to improve the theory of nuclear forces, and others. These hadrons are not sta- 
ble, with for instance A— )■ N + tx or p — > n + n, but were named hadrons as well, baryons 
or mesons. 

The A, at first sight, appears as a consequence of the interaction between -k and N, 
and the p as a resonance of the -kit interaction. Chew and his collaborators generalised 
this scenario, and suggested the concept of "bootstrap" or "nuclear democracy" |36 37 [: 
everything is made of everything, and any hadron is both a building block and the 
result of the interaction of the other hadrons. For instance, A is made of Ntt, Ntctc, etc., 
and, as well iV is made of An, etc. This gives an infinite set of coupled equations, of 
which it was hoped one could extract a finite set as a first tractable approximation to the 
spectrum and to the dynamics. The success was, however, extremely limited. Ball, Scotti 
and Wong, for instance, stressed that describing mesons as resulting of the nucleon- 
antinucleon interaction hardly gives the observed "exchange degeneracy" (named after 
the phenomenology of high-energy scattering), the property that an isoscalar and an 
isovector mesons with the same quantum numbers have very often the same mass |38~| . 
In the baryon sector, the next state after the A with isospin 1 = 3/2 and spin J = 3/2 
was sometimes predicted to have 1 = 5/2 and J = 5/2! 

Third came strangeness. New particles were observed in the 50s and 60s, decaying 
weakly though massive enough to decay to existing particles, and produced by pairs 
with strict rules: K + together with A for instance, but never K~ together A. A new 
quantum number, strangeness S, was introduced to summarise the properties of these 
new particles: strangeness is conserved in production processes by strong interaction, 
and thus A(S = — 1) can be produced in association with K + (S = 1), but not K~ which 
has S = — 1. On the other hand, strangeness is not conserved in the decay by weak 
interaction, as A(S = — 1) — > N + n, or A — > p + e~ + v. The weak interaction of strange 
particle was beautifully linked to that involved in ordinary beta decay. 
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2.3 Generalised isospin 

There is of course the exception of the n meson, with a mass of about 0.14GeV, signifi- 
cantly lighter than the mass of the K, about 0.49 GeV, and the exception of light scalar 
mesons with long-standing questions about their structure. Otherwise, one observes 
that strange particles do not differ much form the non-strange ones. For instance, the 
mass of the A baryon is about 1.1 GeV, just slightly above that of the nucleon at 0.94 GeV, 
and the mass of the K*, 0.89 GeV, is close to that of the vector mesons p and u, about 
0.78 GeV. It was thus tempting to put strange and non-strange particles in multiplets 
generalising isospin. Since isospin is built on the SU(2) group, the minimal extension is 
SU(3). Later, it was renamed SU(3)f, to differentiate it from the SU(3) group associated 
with colour. 

A specific model was proposed by Sakata p9j, with (n, p, A) as the building blocks of 
matter, and mesons as baryon-antibaryon pairs. We already mentioned that this picture 
of mesons is difficult to accommodate with the long-range baryon-baryon interaction 
as given in the Yukawa picture. But the Sakata model had also problems with baryons. 
There are too many baryons with low mass. Take for instance the lowest baryons with 
spin J = 1/2. Besides p, n, and A, there is a triplet de singly-strange (S = —1) baryons 
E°, £~) of mass about 1.3 GeV, and a pair of doubly strange (S = —2) baryons 
) with mass about 1.5 GeV. In the Sakata model, they should belong to higher 
representation, and, in a dynamical picture, contain an additional baryon-antibaryon 
pair. This is hard to believe. 

To take care of this problem, Ne'emann and Gell-Mann suggested to keep the SU(3) 
group as the basic symmetry, but to put the known J = 1/2 baryons in an octet represen- 
tation. This is the famous "eightfold way" pD|. The group SU(3) has eight generators, 



instead of three for SU(2) (I± and I3). Each multiplet can be characterised by the dimen- 
sion of the representation (this generalises the 21+1 multiplicity for SU(2)), and two 
generators that commute, which are usually taken as I3 and strangeness S, or equiv- 
alently the hypercharge defined as Y = b + S, where b is the baryon number. SU(3) is 




Figure 1: Octet of the lowest spin J = 1/2 baryons 



rather good symmetry. Its breaking can be described by simple terms, treated at first 
order, which are proportional to some generators of the group. More elaborate mecha- 
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nisms were proposed for breaking SU(3). One example is the famous Gell-Mann-Okubo 
formula for the baryon masses 

M = M + aY + b{I{I+l)-Y 2 /A) , (6) 

from which one gets 

2(iV+~) = 3A + £, (7) 

where each particle stands for its mass, in surprisingly good agreement with experi- 
ment. 

2.4 The success of the eightfold way 

When the eightfold way was proposed, only 9 baryons with spin J = 3/2 were known, 
with a low mass: four A of charge ranging from —1 to +2, thre S* with strangeness 
S = — 1 and two E* with strangeness S = —2. At the 1962 Rochester Conference held 
in Geneva, Gell-Mann pointed out they would fit very well a decuplet representation of 
SU(3), provided the last member does also exist. He named it f2 , as a kind of ultimate 
achievement, in the biblical sense. The masses of the 9 existing members being observed 
to grow linearly with strangeness, i.e., following the pattern 

S* - A = E* - S* , (8) 

known as the "equal spacing rule of the decuplet", it was tempting to extrapolate to in- 
clude Q~ — E* in the equality, thus predicting the mass of the Q at about 1.67 GeV. Some 
physicists were sceptical about the possibility of producing and detecting easily the Q~. 
It was nevertheless seen in an experiment at Brookhaven led by Nick Samios, with ex- 
actly the mass predicted by Gell-Mann. This was at the end of 1963, and published in 
1964 |4~T| , which otherwise was one of the best vintages ever in Burgundy! 

2.5 The fundamental representation: quarks 

The pseudoscalar mesons (n, 77, K, K) were accommodated in an octet and a singlet. As 
there is presumably mixing between the singlet and the isoscalar member of the octet, 
it became customary to talk about "nonet". The same holds for the vector mesons. See 

We thus had already in the hadronic world octets, singlets and decuplets, but no 
triplet, which corresponds to the fundamental representation! Gell-Mann proposed that 
the fundamental representation is populated by three yet not discovered - or hypothet- 
ical - particles. He was of course fully aware that any representation can be built by 
combining the fundamental representation, 3, and its conjugate, 3, in the same way that 
any value of the spin can be built by adding elementary spins 1/2. 

The word "quark" was taken from the sentence "Three Quarks for Muster Mark" in 
Joyce's Finnegans Wake. In another context (see next subsection), the word "ace" was 
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Figure 2: Decuplet of the lowest spin J = 3/2 baryons. The Q was missing till its 
discovery in 1964. 

suggested, but did not prevail. The individual quarks were sometimes named p, n and 
A, in reference to the Sakata model, but quickly the naming scheme u, d, s was adopted 
by the community. 

The properties of the quarks are summarized in Table [T] and Fig. [4] 

Table 1: Properties of the quarks 
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Figure 3: Octet and singlet of pseudoscalar (left) and vector (right) mesons. 
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The basic reduction of products of representations that are important for the quark 
model are 

3x3x3 = 1 + 8 + 8 + 10, 3x3 = 1 + 8. (9) 

Historians could debate endlessly whether the pioneers considered quarks as a handy 
mathematical tool to build the representation of SU(3), or had already in mind a phys- 
ical interpretation of the quarks as the constituents of the hadrons. Anyhow, this was 
one of the major breakthroughs in physics. 



Y Y 




h 



Figure 4: Triplet of quarks and antitriplet of antiquarks 



2.6 The OZI rule 

Another approach was followed by Zweig. See, e.g., his recollection at the Conference 
Baryon80 g2J or at the celebration of Gell-Mann's 80th birthday (43J. 

The <p meson, of mass about 1.02 GeV, was discovered in 1962 |44l|45j. It is an 



isoscalar, vector meson, like the oj, but with peculiar decay patterns. While it could 
easily decay into pions, it prefers the KK channels, for which the phase-space is mea- 
gre (the K has a mass of about 0.49 GeV). 

Zweig's explanation is that this favoured decay is dictated by 4> meson content. He 
named the constituent "aces", but we will call them quarks to conform with the current 
usage. The idea is that the decay preferentially keeps the existing content. In the modern 
language, the 4> decay is described by the diagrams of Fig. |5j 




7T 



c 



7T 



K 



K 



Figure 5: Decay of the meson: Zweig forbidden (left, centre) or allowed (right) 
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He thereby invented the "Zweig rule", also called Okubo-Iizuka-Zweig (OZI), or 
(A-Z) rule since many authors contributed, from Alexander to Zweig. For a review, see, 
e.g., |46|. As we shall see shortly, a nice surprise was that the rule works even better for 
heavier quarks. 

The rule governing the <p decay was extended to other decays and to reaction mecha- 
nisms. Quark lines should not start from and end in the same hadron, i.e., disconnected 
diagrams are suppressed. Instead quark lines should better link two hadrons in the 
initial or final state. 

The OZI rule then got variants. It is sometimes argued that the dominant processes 
correspond to planar diagrams, while non-planar (but still connected) diagrams are sup- 
pressed. For instance, there are dozens of measured branching ratios for antinucleon- 
nucleon annihilation at rest or in flight. In Fig. |6j the rearrangement diagram (left) 
clearly keeps the initial constituents, while the diagram on the right is more planar. 
Clearly, the former does not produce enough kaons, while the second predicts much 
more kaons than observed. 

Note that antiproton annihilation is an indirect evidence for quarks, but this was 
understood much after the first measurements. If guidance is sought from QED, anni- 
hilation has a rather small cross section, because it is a very short-range process. The 
results obtained at Berkeley for the elastic and annihilation cross-sections of antipro- 
ton scattering on nucleons indicated that the latter is larger. To reproduce these results 
with an empirical (complex) potential, one needs a large size for the annihilation part, 
about 0.8 fm or more. This was a puzzle. One now understands that the proton and 
the antiproton are composite, with a size of the order of 0.5-1 fm. Hence, when they 
overlap, they can rearrange their constituents into quark-antiquark pairs. This is sim- 
ilar to the rearrangement collisions in molecular physics, but has little to do with e + e~ 
annihilation in QED. 




Figure 6: Possible diagrams contributing to baryon-antibaryon annihilation. The hi- 
erarchy of these diagrams requires an extension of the Zweig rule that remains a little 
controversial. 



2.7 First quark models 



Greenberg, in a celebrated paper |47[, tried to understand the structure of baryons made 
of quarks. In this paper, he suggested a kind a harmonic oscillator as a first approx- 
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imation to describe the quark motion. He addressed the problem of the statistics and 
suggested a kind of para-statistics, that eventually became the colour degree of freedom. 

The work of Dalitz was done almost simultaneously. In the summer of 1965, and in 
particular in his lectures at the School of Les Houches fl4"8| , he constructed his first ver- 
sion of the harmonic oscillator model of baryons. As Greenberg, he faced the problem 
of the statistics of the quark. 

Dalitz's work was the starting point of a series of studies about baryons in the har- 
monic oscillator model with contributions by Horgan, Hey, Kelly, Reinders, Gromes, 
Stamatescu, Stancu, Cutosky, etc., culminating with Isgur and Karl. Potential models 
not based on harmonic confinement were proposed somewhat later. The references will 
be given in the section on baryons. 

Note also the contribution by Becchi and Morpurgo about the possibility of describ- 
ing hadrons made of quarks in a non-relativisitic approximation 
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2.8 Heavy quarks 

The physics of kaons has been very stimulating along the years: strangeness led to 
the quarks, the 9 — r puzzle led to parity violation and to K°K° mixing, whose de- 
tailed scrutiny revealed CP violation pV| . Another problem in the weak decay of kaons 



inspired Glashow, Illiopoulos and Maiani |51[, who predicted a fourth quark, named 



"charmed quark" and abbreviated as c, whose mass should not be too high. In some 
processes, diagrams involving a u and diagrams involving a c cancel out. This is the 
GIM mechanism. 

A new spectroscopy was thus predicted, with charmed mesons such as (cu), or 
charmed baryons such as (csu), double-charm baryons, etc. See, e.g., Gaillard, Lee and 



Rosner [52|. Of course, (cc) mesons were predicted as well. 

However, when in November 1974 (this was sometimes called the October revo- 
lution), the J I ip was discovered simultaneously at SLAC and Brookhaven, and the ip' 
shortly after at SLAC, they were not immediately recognised as (cc), i.e., charmonium. 
The surprise was that they are extremely narrow resonances. We now understand that 
the Zweig rule works better and better when the quark become heavier. The spectrum 
of charmonium was completed with various P state, the if)" = ^(3770) which is a D state, 
and after some time, with the spin-singlet states. 

Charmed mesons [53 1 and baryons [54| were discovered as well, and this sector is 



now rather rich, tough the double- and triple-charm baryons still await discovery. 

The charmonium gave a decisive impulse to the quark model in the meson sector. 
Thanks to Regge and others, we had already an idea about sequences of mesons with 
increasing spin J. In the quark model, this corresponds to orbital excitations of the 
quark-antiquark motion. With charmonium, the new feature is that i) the levels are 
better seen, since the lowest states are narrow, ii) there is a clear evidence for the radial 



3 I thank Pr. Morpurgo for a correspondence related to this subject 
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degree of freedom. Explicit quark models were developed to describe the (cc) spectrum. 
They will be reviewed in Sec. |3j 

At the time of the discovery of charm, in the 70s, there were already speculations 
about a symmetry between quarks and leptons. The light quarks (it, d) are the partners 
of (e~, i> e ). The strange and charmed quarks belong in the family of (fj,~, u^). 

Note that the leptons are ahead. The fi~ was discovered in 1936j^]and was completely 
unexpected ("who ordered the muon? asked Rabbi). The r lepton was discovered at 
SLAC in 1977. The partners of the (r~, v T ) pairs were thus anticipated and named (t, b), 
as "top" and "bottom". And a spectroscopy of hadrons containing b or t, or both, was 
predicted. However, at a Conference in Hamburg, a German minister who had some 
knowledge of English, suggested to replace bottom by "beauty". And t was renamed 
"truth quark", but this is not very often used. 

In 1977, Leon Lederman, who missed the charmonium by a small margin, did an 
experiment similar to Ting's, but with a more powerful beam and an improved detector, 
and announced in 1977 the discovery of the T and T' particles, immediately interpreted 
as (bb) bound states 1 55 56\. We shall come back in the next sections about the role of 
flavour independence. Just a word here. Lederman noticed that, within the precision 
of his measurement, T' — T ip' — J/ip, and submitted to local theorists the question 
whether there exists a potential such that changing the reduced mass keep the spacings 
AE unchanged (remember that AE oc m for the Coulomb potential and AE oc m~ 1//2 
for the harmonic oscillator). Quigg and Rosner found that the logarithmic potential has 
this property for all spacings [57|. In fact, the logarithmic potential was already used 
in empirical speculations about the not-yet-discovered bottomonium, but this property 
was not stressed clearly enough |58| . 

Predictions for (tt) were revised continuously, as the mass of the top was pushed 
higher and higher by the negative results of the experimental searches. It was then 
stressed (see, e.g., | |59~| for an early review) that if the top quark is very heavy, it will 
decay before hadronising. 



2.9 Confirmation 

On the way from the early days of the quark model to the recent state of art, there are 
many beautiful and decisive contributions that unfortunately, I cannot review here, due 
to the lack of time. I will just list of few of them: deep inelastic scattering, the parton 
model, the MIT bag model, and its many variants, QCD, its property of asymptotic 
freedom, lattice QCD, QCD sum rules, etc., etc. 



4 In 1947, the discovery of the pion was confused by the decay of pions into muons 
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3 The quark-anti quark model of mesons 



I married them 
Friar Laurence, Romeo and Juliet 



3.1 Introduction 



The quark model of mesons has not been developed "from the bottom up", but mainly 
"from the charm down to light sector, and from the charm to the beauty quark". 

At the end of 1974, when the new particles seen at Brookhaven and SLAC were 
identified as (cc) bound states, explicit models were proposed to calculate the spectrum 



and the radiative transitions |60j - |63j . The potential proposed in [63 1 is know as the 
funnel potential or Cornell potential, and reads, 

V(r) = -- + br + c (10) 
r 

A new era of meson spectroscopy was opened, with explicit calculations of the meson 
properties in models which were first guessed empirically and are now more seriously 
inspired by QCD. 



3.2 Quantum numbers 

The total angular momentum J results from the coupling of the spin S of the quarks 
and their orbital momentum £, from which one also gets the parity (— and charge 
conjugation C = (— l) £+s . The lowest states are given in Tableg 



Table 2: Quantum numbers of the lowest quarkonium states 

25+ Hj l S 3 S X ^ 3 P 3 P X 3 P 2 X D 2 3 D X 3 D 2 3 D 3 
J PC 0-+ 1— 1+- 0++ 1++ 2++ 2-+ 1— 2- 3~ 



Comments are in order 

• some sets of quantum numbers are absent. For instance, a state with J PC = 1 + 
would be exotic, 

• some J PC occur twice. There is the possibility of, e.g., 3 Si — 3 Di mixing, as for the 
deuteron, 

• the above quantum numbers are repeated for the radial excitations, labelled with 
the radial number n. In the literature, one can find either n = 0, 1 ... as in the 
familiar one-dimensional oscillator, or n — £ + 1, 2, . . . as in the Coulomb problem. 
In the latter case, the first P state is labelled 2P, the first D state 3D. We shall adopt 
n = 1, 2, . . ., i.e., a counting where the radial wave function has (n — 1) nodes. 
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3.3 Spin averaged spectrum 

With a potential such as ( |T0"| which is central, without spin dependence, the energy 
depends only on I and n. The radial equation reads 

1(1 + 1) 

— u (r) H — u(r) + m V(r) w(r) = mE u(r) , (11) 

where m = m c is the mass of the charmed quark in the case of quarkonium, and other- 
wise m is twice the reduced mass. For a pure Coulomb interaction, this equation can be 
rescaled to —v"(r) + (£(£ + l)/r 2 — 1/r — e)u(r) = 0, with eigenenergies e = — 1/(4 n 2 ), 
with n = £ + 1,£ + 2, . . .. For a purely linear interaction, another rescaling leads to 
—v"(r) + (£(£ + l)/r 2 + r — e)v(r) = 0, and in the case where I = 0, it reduces to a 
shifted Airy equation: the eigenenergies are the negative of the zeros of the Airy 
function Ai(x), and the n th eigenfunction is just the very same Airy function shifted at 
«W v n (r) oc Ai(r + a n ), the normalised version being 

Ai(r + a n ) 

W = a ■// x • (12) 
Al (0 n ) 

If both the Coulomb and the linear terms are present, one can rescale to a one-parameter 
problem, which can be chosen as 



. A 

-A h r 

r 



Mr) = , (13) 



where A can be expressed in terms of the quark mass m and the strength parameters a 
and b of the Cornell potential. 

If the two-body problem is solved, one can tune the parameters to reproduce the low 
levels of charmonium. This was done by several groups in the late 70s, and the authors 
were able to predict the missing states. 

Note that this exercise did not provide with a sharp determination of the parameters, 
e.g., a, b and c in ( [T0| >, and the mass m c of the charmed quark. Additional constraints 
were used, such as the leptonic widths, and the rates of the 7-transitions. Still some 
flexibility was allowed, and, when the first bottomonium levels were found, the game 
became more challenging: to reproduce simultaneously the (cc) and the (bb) spectra. 

Indeed, even so the interquark potential was not derived from QCD in early quarko- 
nium phenomenology, it was assumed that it is universal, or "flavour independent". In 
QCD the gluons couple to the colour. Hence it is reasonable to assume that the potential 
is flavour independent. 



We already mentioned the logarithmic solution for V(r) |57[. Then the radius scales 
as m _1//2 , and the square of the wave function at the origin, |0(O)| 2 , that enter several 
decay widths, scale as m~ 3 / 2 . As lnr = lim(r a — I) /a as a — > 0, a generalisation consists 
of a power-law interaction V(r) = Ar a + B. Martin used this functional form to fit the 
quarkonium levels, and obtained a rather good fit |64| with a small value of a, i.e., a 
potential close to the logarithmic one. See, also, |65|. 
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A typical charmonium potential is shown in Fig. y\ as well as the reduced radial 
wave function u(r) for the IS and 2S states. In units of GeV for V and GeV 1 for r, it 
reads V(r) = —A/r + 0.2 r — 0.35. It is just for illustration purpose, without any attempt 
to achieve the best fit. 



u(r),V(r) 




Figure 7: Simple central potential — 0.4/r + 0.2 r — 0.35 and two first S-wave levels for a 
quark mass m c = 1.5. Units are GeV for V , GeV -1 for r, and arbitrary for the reduced 
radial functions u(r). 

The corresponding results are listed in Table |3j This should be considered as a kind 
of 0th order starting point which could be refined by tuning the parameters, fitting some 
higher levels, and including some decay properties in the constraints. 



Table 3: Rough fit to the spin-averaged levels of quarkonium (for ID, the experimental 
value corresponds to ip" which is 3 Dx, and the state seen by CLEO (66j and BABAR (67|, 
which is presumably 3 D 2 

bb 

IP ID 2P 
9.87 10.11 10.23 
9.89 10.16 10.26 

In Fig. [8] are shown the levels of bottomonium using the simple potential —0.4/r + 
0.2 r — 0.35 and a mass m& = 4.5 GeV. For the very beginners, its a good exercise to repro- 
duce these numbers. The hierarchy of the excitations corresponds to the observation. 

In Fig. [9j we supplement the previous figure by the charmonium levels, computed 
with the same simple potential. The IS levels are arranged to coincide. The potential is 
tuned to produce about the same spacing for the lowest states. But for the higher states, 
the spectrum becomes dominated by the linear part of the interaction, and the spacing 
is significantly higher for (cc) than for (bb). 



cc 

IS 2S IP ID IS 2S 

Model 3.07 3.68 3.48 3.78 9.47 9.99 

exp. 3.07 3.67 3.52 3.77 9.44 10.01 
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M (GeV) 



10.2 



3S 



2S 



9.8 



IS 



9.4 



2P 
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ID 



Figure 8: Predictions for the bb spectrum with a very simple potential. 



M (GeV) 



10.2 



9.8 



9.4 



3S 
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2P 
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ID 



3.80 



3.40 



3.00 



Figure 9: Comparison of the predictions for the bb (solid lines) and (cc) (dotted lines) 
levels with a very simple potential. 
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3.4 Improvements to the potential 

Perhaps, one should not try to improve too much the simple potential model of mesons. 
Take for instance proton-nucleus scattering in nuclear physics, where the so-called 
Glauber approximation works very well. For years, many corrections have been es- 
timated, and each correction gave a large effect. At the end, one understood that the 
most important corrections cancel each other. There might be something similar in the 
quark model. Anyhow, let us list a few possible improvements. 



3.4.1 More elaborate potentials 

So far, we mentioned the simplest choices for the interquark potential, such as power- 
law, logarithm, or Coulomb-plus-linear. 

The success of these empirical potentials, with a short-range part finite or less sin- 
gular than 1/r, and a confining part less sharp than r can be understood as follows: the 
simplest interaction corresponds to a linear confinement, with a cigar-shaped gluon flux 
linking the two gluons, and a Coulomb-part to which one-gluon exchange contributes. 
But the linear part is smoothed by pair-creation effects; if one tries to increase the sep- 
aration between c and c, a pair of light quark is created, and this results in a softening 



of the interaction [ 68 69 [ . At short-distance, there is the asymptotic freedom, which 
weakens the Coulomb term, see, for instance, the potential by Buchmuller et al. | |70~J or 
Richardson (71). 



3.4.2 Relativistic corrections 

The simplest potential models use the Schrodinger equation. However, the kinetic en- 
ergy turns out not to be very small as compared to the rest mass of the quarks. Models 
have been devised with a relativistic form of the kinetic energy 



Pi 

2m 



— >■ 



\/m 2 + p 2 



m 



(14) 



See, e.g., Basdevant and Boukraa |72| , Godfrey and Isgur |73[, etc. This is of course just 
one step towards a fully relativistic and covariant picture, with retardation effects, etc., 



one of the most serious attempts being the work of the Bonn group [74,75|, which later 
also described the baryons. 



3.4.3 Strong decay of quarkonia 

As already stressed, the miracle of charmonium and bottomonium is the suppression 
of the decay into light hadrons, which proceeds through internal QQ annihilation into 
gluons and subsequent hadronisation of the emitted gluons. Hence the states below the 
relevant (Qq) + (Qq) threshold are narrow, and the states lying above easily decay into 
a pair of flavoured mesons, by a Zweig-allowed mechanism. 
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The most widely used model for describing (QQ) — > (Qq) + (Qq) is based on quark 
pair creation. For, say 

A{q 1 ,q 2 )^B{q 1 ,q 3 )+C{q 3 ,q 2 ) , (15) 

a prescription is given for the creation of the extra pair (q 3 , q 3 ), and the amplitude in- 
cludes the coefficients for the spin-isospin recoupling and the overlap of the wave- 
functions, say 

M = J dr¥(l, 2) 0(3, 3) 3) *(3, 2) , (16) 
where dr is meant for all relative variables. The most popular is the co-called 3 P model 



proposed by Micu |76| and extensively used and developed by the Orsay group |77[ 
and many others. The Cornell model of charmonium |78l is very similar. 



The early applications were mostly for nodeless states, for which (16) gives some 
slight enhancement or suppression of \M\ 2 multiplying pure phase effects. The appli- 
cation to high charmonium states ip^ revealed more dramatic variations. In particular, 
if you consider the three OZI-allowed decays 

each, schematically, call for a certain momentum p in the initial-state wave function. If 
p lies at a node or a bump, the transition is suppressed or enhanced. On the basis of 
simple spin counting, one expects for ( [T7| > that the rates, divided by phase-space, are 
proportional to 

R[DD, D*D + c.c., D*D*] oc 1 : 4 : 7 , (18) 

while the experimental results for ^(4.04) show a clear dominance of D*D*. This led to 
suggest a molecular structure for this ^(4.04) |79| (see the section on multiquarks). In 
fact, the formalism of quark-pair creation shows that the decays into DD or D*D + c.c. 



are suppressed by the node structure of the initial state | [80| . See also, |81 . 82 [. In a 
more detailed analysis, Eichten et al. [83 1 and Fernandez et al. (84j concluded that the 
situation is neither as simple a pure (cc) with nodes or a pure molecule, namely that this 
state is a charmonium with an abundant molecular component^] 

3.4.4 Coupling to decay channel 

The coupling to (Qq) + (Qq) also influence the energy and the internal structure of the 
quarkonia. The quark model is just a first approximation, with the minimal component 
in Fock space, i.e., the first term in 

\ip) = a \cc) + 2^ h i \ c ^iQi) H (19) 

i 

with the second term tentatively saturated by the lowest D^*>D^*> component and the 
D s analogues. In the model by Eichten et al. [78,82], these contributions are generated 



3 I thank F. Fernandez for an interesting correspondence on this point 
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from the main cc component by an explicit operator which creates a pair of light or 
strange quarks out of the vacuum. This provides the states above the threshold with an 
explicit decay width. The branching ratios are found to be in good agreement with the 
observed ones. There is also a dispersive part, i.e., a shift of the energy due to meson 
loops (see Fig. (Tob. If one introduces an explicit mass difference between D and D* , 



these meson loops contribute to the fine and hyperfine splittings. See also [85], and the 



recent contributions by Barnes et al. [86 [, T. Burns |87[, etc., where it is stressed that the 



loops with D, D* , D s , etc., tend to cancel one another. 

£)(*) 





Figure 10: Mesonic loop correction to the charmonium, seen at the hadron level (left) or 
quark level (right) 



However, very close to the DD threshold, there is an effect on the ip' , which can 
couple to DD, while the r]' c only couples to DD* + cc. and D*D*. Hence the — rj c is 
sensitive to the coupling to virtual decay channels |88|[89| . 

We shall come back on these meson-meson configurations in the section on multi- 
quarks. 



3.4.5 Fine structure 

Three of the four members of the IP level of (cc) were discovered rather early, thanks to 
the transitions 

ip' -> Xj + 7 , Xj -> J/ip + 7 , (20) 

with J = 0, 1, 2. The masses have later been measured with a very high precision using 
an antiproton beam |90| . 

In the quark model, they are the spin triplet states 3 P , 3 Pi and 3 P 2 (with a small ad- 
mixture of 3 F 2 ). The formalism for the splittings and for the electromagnetic transitions 
is adapted from nuclear and atomic physics, and is reviewed, e.g., in Jackson's lectures 



atSLACin 1976 |91|. 



The potential is written as 

V(r) + V ss (r) <7i.0- 2 + V ls (r) l.s + V t (r) S 12 , (21) 



where I is the orbital momentum, s = (cri + cr 2 )/2 the total spin of the quarks, and the 
tensor operator is Si 2 = "icn.r cr\.r — ct\.02- If the two last terms in ( [21) are treated at 
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first order, one gets the masses 

M( 3 P ) = M t - 2 (V ls ) - 4 (V t ) , 

M( 3 F 1 ) = M t - (Vi s ) + 2(V t ), (22) 

M( 3 P 2 ) = M 4 + (VJ s )-|(F t ), 

5 

in terms of a an average triplet mass M t and two matrix elements. This can be inverted 
as 

M t = i [M( 3 P ) + 3 M( 3 P0 + 5 M( 3 P 2 )] , 
(V^) = ^ [-2 M( 3 P ) - 3 M( 3 Px) + 5 M( 3 P0] , (23) 

(V t ) = ^ [2M( 3 P ) - 3 M( 3 P0 + M( 3 P!)] . 

Note that if the spin-orbit and tensor terms are treated beyond perturbation theory, 
leading to the computed masses M( 3 Pj), and if M t still denotes the computed result 
with only V(r) + V ss (r), then one gets 

M t >\ [M( 3 P ) + 3 M( 3 Px) + 5 M( 3 P 2 )] , (24) 

i.e., M t is above the naive weighted average. This should be kept in mind, when esti- 
mating the sign of V ss by comparing the spin-triplet states with the spin-singlet state 
governed by V(r) — 3 V ss (r) (see below the paragraph on 1 P 1 ). 

From the date of the Particle Data Group |92| on the IP levels, of charmonium, ( |23] > 
gives 

M( = 3.525, (Vis) = 0.035, and (K)=0.01GeV. (25) 
For the IP level of (bb), the result is 

M t = 9.900, (V u ) = 0.014, and (V t ) = 0.003 GeV , (26) 

and for 2P 

M t = 10.260, (Vis) = 0.009, and (V t ) = 0.002 GeV . (27) 



3.4.6 Hyperfine splittings 

This is a more delicate issue on the experimental side. The simple mechanisms explain- 
ing the fine structure (see next section) also predicted some hyperfine splitting, with a 
pseudoscalar state named r\ c about 0.1 GeV lighter than its vector partner J/ip,. The r\ c 
is the 1 S state of (cc), with n — 1. It was thus anticipated that the Ml transition from 
J/ip to t] c would not have a very large probability [93 94[. 



A first candidate for r] c was claimed in Germany |95[, with a surprising shift S = 
m(J/iJ}) — m(r] c ) ~ 0.3 GeV. Leutwyler and Stern [ 96 [ claimed that such a large shift was 
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an unavoidable consequence of the relativistic character of the dynamics. In contrast, 
most potential builders predicted a smaller 8, once the first r\ c was announced, recog- 



nized to face difficulties to accommodate such a large value of 8 [97 - 101 [. 



The German r\ c was not confirmed 1 102 [, and instead, a more plausible one was seen 
at 2.980 GeV [ 103 104) , and confirmed in other places, including LEP and antiproton- 
proton collisions. The corresponding splitting is 8 ~ 0.1 16 GeV. 

The r] c (2S), sometimes named rj' c , was then predicted with a splitting 8' = ip' — rj' c of 
about 80 MeV. The ratio 8'/ 8 is mainly governed by the ratio R = |0 2 (O) 1 2 /l<M°) 1 2 of the 
wave function at the origin. This ratio can be estimated in potential models (for instance, 
this ratio is exactly R = 1 for a purely linear potential). It can also be deduced from the 
ratio of the leptonic widths. The i] c (2S) was unfortunately forgotten in the elaborate 



Cornell model with explicit account for the coupling to the virtual decay channels |82 [. 



This omission was repaired in [88,89], and the effect of meson loops found to be sizeable 
The r] c (2S) does not couple to the nearby DD, but the ip' does, and is pushed down. This 
reduces 8'. The r] c {2S) was eventually found at BELLE, BABAR and CLEO, in B decay, 



in 77 spectra and double-charmonium production [105 -107[. It means that if one looks 
at e + e~ — > J/ip + X and believes the Zweig rule, X is dominated by (cc), and, indeed, 
already known states have been seen in X as striking peaks. The i] c (2S) was confirmed 
in some other experiments. The splitting is ip' — rj c = 0.049 GeV. 

From the very beginning, it was suspected that the h c , corresponding to the 1 P i state 



of (cc), will be very difficult to produce. See, e.g., Renard [ 108 [. In most models, h c is 
predicted to lie very close to the centre of gravity of the spin- triplet, M t . A first indi- 
cation was provided by the R704 experiment at CERN, the last experiment at the ISR 
ring (109). The h c was later found at CLEO (Cornell) (ITO) and E835 (Fermilab) (111) . 
The mass of h c coincides with the centre of gravity of the triplet, and there is presumably 
a cancellation of several small terms. 

The analogues in the (bb) family were also hard to find. We have now r]b, hb(lP) and 



hb(2P) in the particle listings [92 [. Very recently, CLEO |112[ reported a measurement 
of r/ b (2S) and a confirmation of 7] b (lS). The results are 

5(15) = 0.067, 5(25) = 0.049 h b (lP) = 9.898 , h b (2P) = 10.260 GeV . (28) 

Two remarks are in order: 



The value of 8(1S) for (bb) was predicted by HPQCD (113J to be 70 ± 9 MeV. If 
you cannot afford lattice QCD, you can attempt a poor-man estimate from (cc). 
Assuming a logarithmic potential, with the distances scaling as m -1 / 2 , and a hy- 
perfine splitting proportional to \i(j(0) 2 \ 2 /m 2 , one gets[] 



S(lS) a = 8(1S) 



cc I I 

m c J 



-1/2 



65 MeV 



(29) 



6 I remember doing this computation using the cell phone of Kamal Seth in a fancy restaurant of Mu- 
nich 
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2. The ratio of 6 (IS) to 5(2S) drops from about 2.3 to about 1.4 in (cc) to (66). This 
reflects how anomalously high is the r) c (2S) due to its vicinity of the OZI threshold. 



3.4.7 The origin of the spin-dependent forces 



The phenomenology of the spin-dependent forces in charmonium was greatly inspired 
by QED and by nuclear forces. 

It was early stressed that the one-gluon-exchange contribution to the cc interaction 
is very similar to the Coulomb interaction in QED, and thus is associated to the same 
spin corrections. The name "chromomagnetism" was given to the analogue of the mag- 
netic interaction. We shall come back on the important role of chromomagnetism in the 



systematics of hyperfine splittings in ordinary hadrons [ 114 1 and in speculations about 
possible multi quarks | |115J . 

The one-meson-exchange picture of the nuclear forces, which was very popular in 
the 60s and early 70s. The nucleon-nucleon potential was written V = V^, where each 
Vi corresponds to a single partial wave in the t-channel, and thus gives a very specific 
spin dependence in the s-channel. 



Schnitzer |97 98 [, for instance, analysed the splittings assuming a vector meson ex- 
change, and a scalar confinement. Thus the 1/r terms gives spin-spin, tensor and spin- 
orbit terms similar to the terms describing the fine and hyperfine structure of atoms, 
or similar to the terms associated with u exchange in nuclear forces. The term propor- 
tional to r gives a negative spin-orbit force. Several variants were discussed by various 



authors |100l|116pl7| . 

The method to reach full consistency of this approach was explained by Gromes 
(see, e.g., ]118[ for references) and Eichten and Feinberg [119 [. The problem is that the 
effective potential in the Schrodinger equation contains two types of contributions: an 
intrinsic spin dependence due to the nature of the exchanged object, and a relativis- 
tic effect known as Thomas precession. Even a scalar interaction, when reduced to an 
effective non-relativistic potential, contains a central term and a related spin-orbit one. 



3.4.8 Orbital mixing 

The states with natural parity, except 3 P , contains two partial waves, which can mix. 
Let us consider J PC = 1 . The wave-function reads 



if, 



u ( r ) |3 g j + M r )i3 



Di) 



and after some algebra, it can be shown that the coupled radial equations read 

u"(r) 

ZT^ + V c (r)u(r) + V8V t {r)w{r) = Eu{r) , 

w"(r) 



m 



+ 



m 



+ V c (r)-3V ls (r)-2V t (r) 



w(r) + v8 Vt{r) u(r) = Eu(r) 



(30) 



(31) 
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first derived many years ago for the deuteron. Here Vc = V + V ss is the spin-triplet 



central potential. Solving (31 1 requires a regularisation of V t {r) which behaves as r 3 at 



short distances in most models. Often the result, for e.g., tp" is written as |120[ 

^(3770) = a\ 3 D 1 ,n= 1) +6 1 | 3 S 1 ,n = 1) +6 2 | 3 Si,n = 2) + • • • 



(32) 



where the states in the r.h.s. are obtained by neglecting the coupling of the two equa- 
tions. One could debate whether b 2 is larger than b x because tp' lies very close to tp", or 
b\ is larger due to the favourable node structure of J ftp and tp". This is not obvious: in 



atomic physics, the analogue of (32 ) requires a summation over all bi (in the discrete and 
continuous spectrum). See, e.g., |121| 122[. 

The orbital mixing is responsible for the coupling of ^(3770) to e + e~ that made its 
observation possible. It might also influence the pattern of its hadronic decays. See, e.g., 
the discussion by Rosner [ |120[ . 

Very likely, a good fraction of the orbital mixing is due to the coupling to decay chan- 



nels. For instance, the coefficients in ( |32| > are estimated in |89| , which uses an improved 
version of the Cornell model. 



3.5 Summary for heavy quarkonia 

In Fig. [TTJ reproduced from pl| with the kind permission of the author, are shown the 
levels of charmonium. The states labelled X, Y, Z will be discussed in another section. 
It is remarkable to see exactly the levels predicted by the quark model. On the left, the 
spin singlet states, with relative angular momentum i = for the two r] c , and i = 1 for 
h c . In the middle, the various states with spin triplet, and I — 0, 1, and even one example 
with i = 2, and the various recoupling of spin and orbital momentum. The hierarchy 
of orbital vs. radial excitations, and the patterns of fine and hyperfine splittings are 
now well described with simple potentials that are now better understood from the 
underlying QCD. 
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Figure 11: Levels of charmonium, and transitions among them (from |31 [). 
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In Fig. 12 we reproduce the lower part of the charmonium spectrum, and the bot- 
tomonium spectrum, borrowed from the Review of Particle Properties [92 J. 

The spectra are very similar. This is due to the property of flavour independence. 
The differences are also understood, due to explicit mass-dependence of the fine and 
hyperfine terms. 



3.6 Light mesons 

The simple quark model is in principle not applicable to light quarks, but some reck- 
less attempts were encouraging. For instance, Martin | |64[ managed to fit (bb), (cc), (ss 



and (cs) with a single potential. The game was extended further by Bhaduri et al. [ 123 [, 
however their uniform regularisation of the hyperfine interaction gave poor predictions 
for J/ift — r} c and similar splittings. Their model was refined by several authors, see, 



e.g., [124 [, with a regularisation of spin forces that depends on the system under con- 
sideration. 

Isgur et al. |1251 have explained why there are some good reasons to link potential 



models, with a minimal amount of relativity, to the underlying QCD. 

Let us give two examples. The first positive-parity excitations of (qq) with 1=1 
are a (980), &i(1235), ai(1260) and a 2 (1320). In the quark model, they correspond to the 
partial wave 3 Po, 1 Pi, 3 Pi and 3 P2- The level order of this IP multiplet, and even the 
pattern of spacings is very similar to what is observed in the charmonium IP states. 

If one looks at meson with high spin J and plot the square mass M 2 against J, one 
finds an almost perfect linear behaviour for this "Regge trajectory". If a linear poten- 
tial is plugged into a Schrodinger equatiion with relativistic kinematics, this linear be- 
haviour is reproduced. 



3.7 Heavy-light mesons 

This is an even more dangerous playing ground. Remember, for instance, that an elec- 
tron is more relativistic in the hydrogen atom, (p, e~), than in the positronium atom, 
(e + , e - ). There is nevertheless a very interesting spectroscopy of D and B mesons. 

In first approximation, the reduced mass governing the internal dynamics is the 
same for (cq) and (bq) mesons, so the excitation spectra should be very similar. This is 
rather obvious in potential models, but it remains valid much beyond this framework. 



This is one of aspects of "heavy quark symmetry" [126- 128 [. 



3.8 Mathematical developments 

The quark model motivated some studies about the properties of the Schrodinger oper- 



ator with confining interactions. See, e.g., [25 129 [ for a review and references. 



Among the subjects, let us mention the following. 
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3.8.1 Level order 

A sufficient condition has been derived that transforms the Coulomb degeneracy 

E(l,£) = E(2,£-l) = ■■■ = E(£ + 1,0) , (33) 
into a series of inequalities: the condition is related to the sign of the Laplacian of the 



potential. As A V > for the funnel potential < pT3] > or any plausible interquark potential, 
this theorem explains the observed IP < 2S ordering in quarkonium, and 2P < 3S in 
bottomonium. Note that AV < for the last electron of an alkalin atom and AV > 
for a muon inside a nucleus: this explains the observed breaking of the Coulomb 
degeneracy for alkalin atoms and muonic atoms. 

Similarly, the sign of the second derivative of the potential as a function of r 2 governs 
the breaking of the HO degeneracy, 

E(l,£) = E(2,£-2) = E(2,£-4) = ■■■ . (34) 

This explains why 2S < ID, i.e., m(^p') < m(ip"). 

There are also results on the spacing of levels as a function of the radial number. 
For the harmonic oscillator, there is a remarkable equal spacing, E(n + 1, £) — E(n,£) = 
E(n, £) — E(n — 1, £). For a linear potential, the spacing decreases. 



3.8.2 The wave function at the origin 

Many decay widths are proportional to the square wave function at the origin, which 
for S-states, reads 

Pn=|$n(0)| 2 = -^<(0) 2 . (35) 
4 7T 

if the reduced radial wave function is arranged to be real. The studies on quarkonium 
gave the opportunity to remember the rule by Schwinger (see, e;g., [ 129J ) 

r°° c\v 

u'(0) 2 = 2fi / ^-u 2 (r)dr. (36) 
Jo dr 

This explains why the p n are independent of n for a purely linear interaction, and, con- 
versely, why the normalisation factor in ( 12 ) is given by the derivative. 

There are sufficient conditions ensuring that p n increases or decreases as a function 
of n (25). 



3.8.3 The dependence upon the masses 

If flavour independence is accepted as a strict rule, one has to study how the energy 
levels and the wave functions evolve when the reduced mass is modified. 

As the kinetic energy term is p 2 /(2 jj), obviously E \ as // /•. The effect is pro- 
nounced when going from (cc) to (bb). On the other hand, the reduced mass governing 
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(Qq) does not change much from Q = c to Q = b. Hence the OZI threshold (Qq) + (Qq) 
has almost the same energy for Q = b and Q = c. This explains why the T(1S) is more 
deeply bound than J/tp with respect to this threshold, and why there are more narrow 
states in the bottomonium spectrum than in the charmonium one. 

One can go further and notice that in a strictly flavour-independent interaction, the 
inverse reduced mass enters linearly There is a known result that for an Hamiltonian 
depending linearly on a parameter, say H = A + XB, the ground state energy e , or 
the sum of first energies ^2f =0 e% is a concave (or convex upward) function of A. Once 
the constituent masses are added, this implies for the ground state of each angular- 
momentum sector, 

(QQ) + (qq) < 2(Qq) , (37) 



an inequality independently derived by Martin and Bertlmann [25 130 1 and Nussinov 
|131[ 132| for potential models, and extended much beyond this framework [ 132[ 133 [. 
Early applications introduced dangerously light quarks, but still, the inequality is satis- 
fied for the spin-triplet states of (cc), (ss) and (cs). Safer is the prediction that the (spin 
averaged) ground state of (be) is above the middle of (cc) and (bb). 

Note that in atomic physics, it is reasonable to believe that the Born-Oppenheimer 
potential that binds di-atomic molecules remains unchanged if the nuclei are replaced 



by isotopes. Thus ( [37] ) holds for H 2 , D 2 and HD, and this governs their abundance at 



thermodynamic equilibrium [134[. 



3.8.4 Concavity with respect to the spin coefficients 

In the same vein, one can explain the sign of the error when treating the spin corrections 
to first order. For the hyperfine effect of S-state alone the singlet and triplet Hamiltoni- 
ans are respectively, in an obvious notation, 

H s = H — 3 V ss , H t = H + V ss , (38) 

so that the fictitious eigenstate of H alone lies at 

E > \e s + -Ek . (39) 

Note that the inclusion of the tensor force would slightly lower E t/ with 3 Si becoming 

:! S: 

It is straightforward to extend this to the case of two parameters. Consider the 
generic spin-triplet Hamiltonian 

Hj = H t + ajV ls + (3jV t , (40) 

so that H + 3 Hi + 5 H 2 = 9 H t , then the fictitious spin-triplet state free of spin-orbit and 
tensor corrections lies at 



E t > [E J=0 + 3 E J=1 + 5 E J=2 \ /9 , 



(41) 
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i.e., above the naive centre of gravity, as discussed earlier in this section. Thus if the x Pi 
is found at the location of the naive of gravity, it has received some attractive spin-spin 
downward shiftfl 



4 Baryons 

Tires faciunt collegiun^ 

. „ T , , , . Latin sentence 

4.1 Introduction 

The quark model of baryons was developed in the 60s, in particular by Dalitz and his 
collaborators, and several other groups, at a time when many nucleon and hyperon 
resonances were already known. 

These early studies, and the more recent work by Isgur and Karl, Gromes et al., 
Stancu et al., Cutosky et al., cannot be separated from the harmonic oscillator (HO) 
model, which provides a powerful classification scheme, and an efficient tool to imple- 
ment the symmetry constraints. This is reviewed, e.g., by Hey and Kelly | |135[ |. 

The HO model has reached a high degree of refinement, with all the possible correc- 
tions carefully listed, and their tentative interpretation in QCD (one-gluon-exchange). 

There has also been attempts to solve the three-body problem, using known tech- 
niques developed in atomic or nuclear physics, such as Faddeev equations, hyperspher- 
ical harmonics or variational methods. Then any potential can be used, for instance a 
funnel interaction with a superposition of Coulomb and linear contributions. 

Note that in most early studies, it was implicitly assumed that the potential is pair- 
wise, and the question has been addressed of the link between the interquark poten- 
tial in baryons and the quark-antiquark potential in mesons. This will be discussed in 
Sec. 4.12[ In fact, it has been anticipated for many years, that the confining interaction 



is of three-body nature, i.e., depends simultaneously upon all relative distances. This 
makes the calculations technically more involved, but makes a link with the lattice QCD 
which favours such potentials [ 136J . 



The role of relativistic effects is of course important. In the latest works by Isgur 
et al., a relativistic kinematics was adopted | |137[ . There are even models based on the 



Bethe-Salpeter equation [138]. 



4.2 Jacobi coordinates 

In the case of equal masses, a set of Jacobi coordinates that diagonalise the kinetic energy 
is (see Fig. [13]) 

„ ri + r 2 + r 3 2r 3 ~r 1 -r 2 

R = g , x = r 2 - ri , y = -j= , (42) 



7 I remember discussions on this point at a Workshop organised in Genoa, to celebrate years of charm 
physics with antiprotons, and later with Kamal Seth. 



8 Three makes a company 
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Figure 13: Jacobi coordinates for the relative motion of the quarks inside a baryon 



where the factor in y makes it easier to implement the symmetry constraints, as ex- 
plained shortly. Then, once the centre of mass motion is removed, the intrinsic Hamil- 
tonian reads 

h = ^ + ^ + V(x,y), (43) 
m m 

where the potential energy, already written as translation-invariant, has to be scalar. 

In the case of two different masses, say (m, m, M), R is modified, but one can keep 
x and y as above. Then the intrinsic Hamiltonian becomes 

2 2 

h = ^ + ^ + V(x,y), (44) 

with the reduced masses now given by 

[i x = m, /i" 1 = (m" 1 + 2M' 1 )/3. (45) 

If the three masses are unequal, the Jacobi coordinate y becomes y oc (mi + m 2 )r 3 — 
m! t*! — m 2 r 2/ and it is straightforward to derive the reduced masses [i x and [i y as func- 
tions of the individual masses m^. 

Later in this section, we shall elaborate more on the potential energy V. Let us just 
mention here two popular choices. The harmonic oscillator corresponds to V = K{r\ 2 + 
r 23+ r ii)/ where = jr^l and r ^ = Tj — Vi. It becomes V = 3 K{x 2 + y 2 ) /2 for (m, m, m) 
and (m, m, M), and a more complicated quadratic form of x and y for (m 1; m 2 , m 3 ). The 
pairwise models reads v (r 12 ) + f (r 23 ) + f (r 31 ). 

4.3 Permutation symmetry 

Schematically, the wave function of a baryon in the quark model is 

V = <iP(x,y)ip s Tpiip c , (46) 



with orbital, spin, isospin and colour parts, and this wave function has to be antisym- 
metric with respect to the permutation of the identical quarks. 
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If there are two identical quarks, as in H = (ssd) or in A = (uds) in the limit where 



isospin symmetry is exact, then each term in (46) is either symmetric or antisymmetric. 
For the ground state of A, with isospin 1 = 0, the colour, spin and isospin parts are 
antisymmetric and the orbital part symmetric. This means spin for the light quarks, 
and thus s = 1/2 for the three quarks, isospin 1 = 0, and ip(x,y) even in x. For instance, 
in the harmonic oscillator, ip(x, y) oc exp[— ax 2 — by 2 ]. 

There are two type of P-states: one where the degree of freedom associated to A is 
excited, with the same symmetry pattern as the ground state. Another wave where the 
light quarks have their spin coupled to s qq = 1, and the orbital wave function ip(x, y) is 
odd in x. In the harmonic oscillator, the orbital wave functions are of the type 

y exp[— a x 2 — by 2 ] , x exp[— a x 2 — b y 2 ] . (47) 

For the former, the quark spin s = 1/2 and the angular momentum £ y gives a baryon 
spin J=l/2orJ = 3/2. For the latter, both s = 1/2 and s = 3/2 are possible, and thus 
a variety of values for the spin J. Usually, the spin orbit splittings are small in baryons 
(see below), and the states with various combinations of spins and angular momenta are 
nearly degenerate. One exception is the A(1405) — A(1520) pair, which has motivated 
an abundant literature. The most plausible explanation is the coupling of one of these 
states to a nearby threshold. 

For three identical quarks, imposing the constraints of antisymmetrisation is more 
delicate. The colour coupling 3x3x3— > 1 is antisymmetric, thus the rest of the wave 
function should be symmetric. This was, indeed, one of the motivations for introducing 
colour. For the Cl~ , this is rather simple: the spin wave-function is symmetric, isospin 
does not exists, and the orbital wave function is symmetric. For instance, xp(x,y) = 
exp[— a(x 2 + y 2 )] in the HO case. This is similar for the A(1232) multiplet, as the isospin 
coupling 1/2 x 1/2 x 1/2 — > 3/2 is symmetric. 

For the first parity excitations of with J p = (l/2)~ or (3/2)~, the spin and the 
orbital wave functions are neither symmetric nor antisymmetric, but their combination 
is symmetric. One cannot escape here to introduce the concept of mixed symmetry. The 



prototype is given by the Jacobi coordinates x and y of ( |42| ). The permutation P 12 reveals 
x odd and y even, but a circular permutation gives 

P^[y + ix]=j[y + ix], (48) 

where j = exp(2z7r/3), as usual. Obviously, if two pairs, say z = v + iu and Z = 
V + iU, have the same permutation properties as z = y + i x, then their coupling give 
a symmetric term, and antisymmetric one, and a new pair of mixed symmetry (again in 
the complex notation), respectively 

^e[zZ*] = uU + vV , 

Sm[z Z*] = v U - u V , (49) 
[z Z]* = (u U - v V) - i (u V + v U) . 
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In the case of the P-states of we have a pair of mixed-symmetry orbital wave 
functions, {ip x ,ip y }, with either x or y excited, which for the harmonic oscillator read 
{i^xi ^y} — {x, y} exp[— a(x 2 + y 2 )], and a pair of spin 1/2 wave functions 

s x = [tit - itt] , s v = ^ i 2 tti - tit - itt] , (so) 



which bears an obvious analogy with the Jacobi variables (42 ). With a spin-independent 
interaction, the P-states of have a wave function of the type [i/j x S x + i/j y S y }/\/2. 

For the nucleon ground state, the orbital part is symmetric, but the spin and isospin 
parts are both of mixed symmetry, with ( [501 for the spin, its analogues {I x , I y } for 
isospin, and the combination [S x I x + S y I y \Jy2 being symmetric. 

The most interesting example of baryon in this respect consists of the antisymmet- 
ric spin-isospin wave function [S x I y — S y I x ] / a/2 which requires a fully antisymmetric 
orbital wave function. This wave function is excited in any interquark separation. For 
instance, in the specific case of the harmonic oscillator, it reads x x y exp[— a(x 2 + y 2 )], 
i.e., £ p = 1 + . This is one of the states predicted in the three-quark model and absent in 
the quark-di quark model. See below the section on diquarks. 



4.4 Solving the three-body problem for baryons 
4.4.1 The harmonic oscillator 

After rescaling, the linear oscillator reads h = p 2 + x 2 with energies e n = 1 + 2 n and 
wave functions (po(x) oc exp(— x 2 /2) and (p n (x) oc [— d/dx + x] n (po(x). 

For equal masses m, the three-body problem with V = 2 K/3(rf 2 + rf 3 + r^) is gov- 
erned by the Hamiltonian 

2 2 

H(m,m,m) = ^ + ^ + K(x 2 + y 2 ) , (51) 
m m 

which corresponds to a sum of two independent spatial oscillator, or six linear oscilla- 
tors. The energies E = (6 + 4 n x + 2£ x + 4 n y + 2 i y ) ^jK/m have a cumulated number 
of excitations iV = 2 n x + £ x + 2 n y + £ y . Here, the radial and orbital excitations in each 
variable are counted separately. 

There is a considerable degeneracy of the levels of given N > 2, which contains con- 
tains several J p states. The literature on baryons within HO has developed a system- 
atic notation which is used beyond this model. In particular, the multiplets are denoted 
[d, £ p }( n \ where d is the spin-flavour multiplicity, £ the total orbital momentum, and (n) 
the radial occurrence (n = is omitted). For instance, the ground state is labelled as 
[56, + ], and the Roper resonance and its partners [56, 0]'. The ground state contains the 
octet of Fig.[l]and the decuplet of Fig.|2j i.e., 2x8 + 4x10 = 56 states. The first orbital 
excitation is [70, 1~], and the first state with a full antisymmetric orbital wave function 
is [20, 1+]. 
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The case of two different masses reads 



2 — 

H(m,m,M) = ^ + I ^ + K(x 2 + y 2 ) , (52) 
m fx 



with [i given by (45). There is still an exact decoupling, and the energies are 



E{m, m, M) = J— (3 + 4 n x + 2 4) + J — (3 + 4 n y + 2 i y ) . (53) 

V tTl y fl 

Here the lowest excitation is explicitly associated with the heavier quark(s),Le., this is 
an excitation in the v ecx variable. This point will be further discussed in the section on 
double-charm baryons. 

The case of three different masses is also solvable exactly, though this is slightly less 
known J^] Once the kinetic energy is written as (44), one can use a rescaling x — > x/y/JI^ 
and y — > yj J~iiy. The Hamiltonian becomes 

H(m 1 ,m 2 ,m 3 ) = p 2 x + p 2 y + Ax 2 + B y 2 + 2C xy , (54) 

and one is left with the diagonalisation of a 2 x 2 matrix. 

4.4.2 First-order perturbation around the harmonic oscillator 

We now turn to more general models for the interaction. The most spread method 
consists of writing the potential as 

V = K(x 2 + y 2 ) + 5V , (55) 

and treat the second term as a perturbation. If K is optimised, this is the first step 
of a variational expansion, to be discussed below. This is of course not very accurate, 
especially at short distances. If 5V tentatively pushes down the Roper resonance (to be 
discussed shortly) from the N = 2 level of the harmonic oscillator to the vicinity of the 
N — 1 level, then perturbation theory is in principle not applicable, as it requires that 
the shifts are small compared to the initial spacings. 

Isgur and Karl ]139[ have noticed the nice patterns of the N = 2 multiplet splitting 



when perturbed by a pairwise SV. It is illustrated in Fig. 14 Gromes and Stamatescu 
]140[ pointed out that the upper part survives if the HO is perturbed by a symmetric 
3-body potential, while the [56, + ]' decouples. The same conclusion was reached by 



Taxil et al. 1 141 [, in the more general situation of a nearly hypercentral potential. 



The upper part of Fig. 14 fits rather well the spin-averaged values of the observa- 
tions. 

This pattern and the one for the N = 3 multiplet have been further discussed by 
Taxil et al. | [i"4"Tfl , Stancu et al. Jl42l , etc. 



9 I remember discussions on this point with J.-L. Basdevant, and with D.O. Riska 
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0.2A 

0.2A 
0.1A 



0.5A 



[20, 1+ 

[70, 1+ 

56,2+ 
70,0+ 



[56, 0+]' 



Figure 14: Splitting pattern of the N = 2 multiplet for a two-body perturbation of 
the harmonic oscillator. [Check whether the 2-body character is required.] There is 
an overall shift (not shown), and, in addition, a shift 0.1 A, 0.2A, 0.3A or ±0.5A of the 
sublevels. 



4.4.3 Converged variational methods 

Let us consider the ground state with equal masses, but the generalisation to unequal 
masses is straightforward. The above method can be read as the first term in a system- 
atic expansion 

■n 

where n denotes the number of excitations and other quantum numbers associated to 
the solutions 4> n (x, y) of the harmonic oscillator. The matrix elements can be calculated 
with subtle recursion relations. The convergence as a function of the maximal n is ex- 
cellent for the energy but elusive for the short-range correlations. 

A variant is the so-called Gaussian expansion, widely used in theoretical chemistry 
and other quantum problems. It reads 

*(x,y) = J2^ exp[-(a^ 2 + 6 l? / 2 + 2c t x.y)] . (57) 

i 

Note that the permutation symmetry is violated in an individual term if cij ^ hi and 
Q 7^ 0, but is restored in the summation either explicitly (correlated Gaussians) or by the 



numerical adjustment. Different techniques are used [143 144] to avoid redundancies 



in the expansion and guide towards an efficient optimisation. 



4.4.4 Faddeev equations 

The Faddeev equations for the three-body problem have been written for short-range 
potentials in momentum space. Later, they have been translated into position space 
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and adapted to the confining interactions. They imply a summation over the internal 
angular momentum, but even with the crudest truncation, they provide an excellent 
wave function. The first paper (to my knowledge) where the Faddeev equations are 
applied to baryons is ]145[. See, also, 1 146 [. 



4.4.5 Hyperspherical expansion 

The 6-dimensional vector {x, y}, car be written (g, f2 5 ) in spherical coordinates. The five 
angles fi 5 include the physical angles x and y and arctan(y/x). The method consists of 
expanding ^(x, y) into partial waves in this space. For a harmonic oscillator, V oc g 2 , 
and one recovers the familiar results. Otherwise, one gets an infinite set of coupled 
equations 

- u'[ L] {g) + ( L + 3 /2)(£ + 5/2) um + J2 V [LUL y(g) u [L y(g) = Eu [L] (g) , (58) 

6 [L]' 

where [L] denotes the "grand" angular momentum and its associated magnetic num- 
bers. 

Keeping only one term in the expansion, the lowest hyperspherical harmonics, gives 
the hyperscalar approximation, which is already very good, as often stressed by the 
Genoa group [ 147J . The first application of this method to baryons is by Guimaraes et 



al. 1 148 [, and Hasenfratz et al. 1 149 1. The convergence was studied in the thesis work by 



Taxil, and is summarised in |146| |1501. 



The splitting pattern of Fig. 14 has been re-analysed, with the starting point being 
now the hyperscalar approximation. The upper part is not modified, but the radial 
excitation now decouples. This gives more freedom to push this state down, but this is 
not sufficient. Hogaasen and I have noticed that the same hyperspherical potential V 00 
governs the ground state, its radial excitations and the first negative parity excitation. 
The effective angular momentum in the radial equation is £ = 3/2 for the former, and 
£ = 5/2 for the latter. Then for any reasonable model for the interquark interaction, the 
radial excitation is above the orbital one ]151[ . The result relies on the local character 
of the interaction and remains valid for local three-body terms. It does not hold if the 
potential is velocity or spin dependent^ 



4.5 Light baryons 

This is the first sector to which the quark model was applied. One cannot compare 
the refinement of the model and the quality of the fit with the analogues for heavy 
quarkonia. Nevertheless, the quark model explains the rough features of the observed 
spectrum, in the sense that any low-lying resonance corresponds to a level predicted in 
the quark model. 

10 I thank W. Plessas for an enlightening discussion on this point, many years ago. 
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Among the problem is that of the Roper resonance. This denotes the state iV*(1440) 
with the same 1 = 1/2 and J p = (l/2) + as the nucleon. Experimentally it lies at about 
the same mass as the orbital excitations of negative parity, perhaps slightly below. In the 
simplest quark model, it is a radial excitation. In the harmonic oscillator, it belongs to 
the N = 2 multiplet and the first orbital excitation to iV = 1. With a standard Coulomb 
+ linear interaction, with an accurate treatment of the three body problem, there is no 
way to make these states degenerate. This is a serious limitation for the naive quark 
model, which misses the relativistic effects and chiral dynamics. For instance a model 
with Goldstone-boson exchange, proposed by Glozman and Riska (T52|, and further 
developed by the Graz group 1 153| 154[ , gets the Roper resonance right. See, also, [ 155 [, 
and the comments by Isgur on this approach [ 156 1. In recent lattice simulations with 
large pion mass (this facilitates the computation, see next lectures), the same ordering 
as potential models is obtained. When the pion mass decreases (and the cost of the 
computations increases), there is a dramatic change in the level ordering, and the exper- 
imental one is recovered (157). 



4.6 Missing states and the diquark alternative 

Conversely, there are several instances where a state predicted by the three-quark model 
is missing in the experimental spectrum. As often stressed by Isgur and Karl (their 
papers in the bibliography and private communications), most experimental data deal 
with pion or photon scattering off a nucleon. This favours states with one quark excited 
and the two remaining quarks nearly untouched, and suppresses states with double 
internal excitation. However, experiments with higher statistics should be able to see 
the missing states. 

A drastic alternative is the diquark models. I should better say the diquark models, 
as there exists different points of view, and also, the model is regularly reinvented by 
authors who do not quote the pioneering papers. For a rather comprehensive review, 
see |I58| . 

Diquarks have been used for many purposes, in particular multiparticle production. 
In spectroscopy, the quark-antiquark model retains only those baryons consisting of a 
compact diquark and a quark. The question is whether the diquark is postulated as an 
elementary constituent or somehow generated dynamically. 



Martin has given a partial answer [ 159 [: for three-quark systems with high orbital 



momentum I, the ground state consists of two quarks mainly in S-wave with the third 
quark carrying most of the orbital momentum, while for low i, the three quarks share 
t symmetrically. This explained, at last, a long standing issue: the Regge trajectories of 
mesons and baryons are linear and have the same slope. 

A digression here about the Regge trajectories. They were invented to explains some 
properties of potential scattering, and later, in a different form, for the high-energy 
scattering of hadrons. In this later context, a Regge trajectory is the squared mass of 
hadrons, as a function of the spin J. The observed trajectories are beautifully linear. 
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For mesons, this was interpreted a flux of constant section linking the quark to the anti- 
quark. 

An intriguing observation is that the baryons also belong to linear trajectories, with 
the same slope. This is of course natural in the quark-diquark model, as the energy is 
that of colour 3 linked to a colour 3. We now understand that this occurs spontaneously. 

A question that is sometimes eluded is whether diquark are just an effective degree 
of freedom that simplifies the dynamics of baryons, or a constituent to be used on the 
same footing as quarks in the "lego" puzzle of hadrons. Recently, Maiani and his col- 



laborators described the new X, Y, Z resonances as diquark-antidiquark states 1 160 [, as 
done earlier for baryonium. Frederikkson et al. | |161| raised the important question for 
the three-diquark systems (though there is a slight confusion with the quantum num- 
bers, the paper remains relevant): do we have a "demon" deuteron made of three di- 
quark, in addition to the ordinary deuteron made of two baryons? 

4.7 Baryons with a single heavy flavour 

Some years ago, this sector was marginal for quark dynamics, with only the ground 
state of (Quu), (Qud), . . . , (Qss). Of course, the measurement of the lifetimes and semi- 
leptonic branching ratios were very exciting, to compare the role of VT-emission, W- 
exchange, W -formation in various flavour configurations. But this is out of the scope of 
this review. 

In recent years, a lot of data came from CLEO, B-factories, Tevatron, etc., and we 
expect many results from LHCb. It has even been suggested that we can learn more on 
light quark dynamics for the {Qqq) systems than from (qqq). This is a domain where the 
data on charm are or will be confirmed by the data on beauty. 

Just to give an example, take f2&. When its mass was published by DO [ 162], it was 



a surprise that — A& was found significantly higher than fi c — A c . See, e.g., [28 



163 [. 



The CDF collaboration got a lower mass |164[, as this is now confirmed by LHCb [ 165 [. 
The mass of this state was predicted by Lipkin et al. |l66 j^j] on the basis of the chro 



momagnetic interaction [ 114[ 115[. See, also, ]167[ 168[. The lesson is that the spectrum 



has regularities, due in particular to flavour independence, and any departure should 
receive an explanation or disappear. 



The spectrum of single-charm and single-beauty baryons are shown in Fig. 15 It 



includes the very recent discoveries at LHC [ 169||170[ . The spectra are arranged so that 



A c and A& coincide. Clearly the excitation energies are very similar. 
4.8 Baryons with two heavy quarks 

It becomes intriguing or even annoying that no (QQq) has been firmly identified yet. 
We have just the results by the SELEX collaboration at Fermilab [ 171[172| , and negative 



results in a few other attempts. As Brodsky often insists in Conference talsk or private 



ii 



I thank Marek Karliner for a correspondence on this point. 
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Figure 15: Comparison of (cqq) and (6gg) experimental spectra. Note that some of 
quantum numbers have not been established, and are the quark model predictions 
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communications, see, e.g., |173|, it is important to detect double-charm baryons. We 
also advocate that one should look for double charm baryons, or even, in the same 
experiments, at "ordinary" double-charm baryons and possible double-charm mesons 
to be discussed in the next chapter. The production of two cc pairs in e + e~ collisions 
occurs at a rate that is somewhat larger than expected. It has been used to analyse the 
mass spectrum of e + e~ — > J/ip + X, where X is mainly cc, thanks to the OZI rule. This 
led to the discovery of r] c (2S). One should reasonably expect that sometimes the two 
pairs hadronise as (ccq) + (cqq) + (cq) + ■ ■ ■ . 

The dynamics of (QQq) baryons is rather rich, with, in the same hadron, the rel- 
ativistic motion of a light quark around a coloured source, as in D mesons, and the 
slow motion of two heavy quarks, as in charmonium. The (QQq) baryons were listed 
in the pioneer paper by Gaillard et al. on charm physics fl52| . The first detailed study 



in potential models was given in 1 174 1 and has been improved by several authors, see, 
e.g., (175) 

At first sight, a quark-diquark approximation seems tempting. Indeed, the average 
distance between the two heavy quarks, tqq, is smaller than the distance r Qq between 
one of the heavy and the light one. However, if the first excitations occurs within the QQ 
subsystem, and in the quark-diquark picture, one should reconstruct a new diquark for 
each of them. 

An alternative approach is the Born-Oppenheimer approximation, as done for H 2 + 
in atomic physics. For a given x, one estimates the light-quark energy and deduce the 
effective QQ interaction^ If this procedure is applied to a simple quark model, one can 
compare the results to the exact numerical solution of the three body problem [174], and 
it turns out that the Born-Oppenheimer approximation is extremely accurate^} Now, 
the effective QQ potential of (QQq) can be calculated on the lattice. 

The problem is therefore where are the first excitations of this system. The question 



has been raised in several recent papers, for instance Cohen, Roberts, etc. [1 76 177\ . The 
answer governs the strategy to be adopted for these systems. 

If the (QQ) subsystem enters the regime where the Coulomb interaction oc a s (M)/r 
dominates, then the spacing will behave as a 2 s M (not a s M as suggested by a misprint 



in [ 177 1). Hence the diquark is frozen in its ground state, and the dynamics is that of a 



light antiquark around a colour 3 source. This is ve ry sim ilar to the physics of flavoured 



mesons (Qq), and this symmetry is promising [176 177 [. 

However, the Coulombic limit of heavy-quark spectroscopy is rather elusive, as 
noted years ago by ]178 [. For instance, if one uses the simple quarkonium potential 



( JlQ[ > and study the level spacing (2S)-(1S) of (QQ) as a function of the quark mass M, 
one finds that it decreases for M < 6 and increases for M > 6GeV. So, for the dou- 
ble charm baryons, the first levels are probably excitations in the diquarks, and can be 
described in the Born-Oppenheimer limit with an effective QQ interaction. 



12 One actually works with the Jacobi variables and thus the recoil effects are properly included 
13 With an energy slightly below the exact one, unlike a variational approximation 
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4.9 Triple heavy flavour 



A classic introduction is a paper by Bjorken |T79| who named this sector "the ultimate 
deal of baryon spectroscopy". This is the analogue of charmonium for baryons, i.e., 
a sector where the heavy-quark dynamics is probed without any complication due to 
light quarks. There are already several predictions for the lowest states. Of course, (bbb), 
(bbc), etc., states are expected as well. 

One of the issues is whether (ccc) has s normal hierarchy of excitations, with first 
the orbital excitation with negative parity, and next the radial excitation which is the 
analogue of the Roper resonance. For light baryons, an anomalous ordering is observed, 
probably due to the chiral dynamics of light quarks, which cannot be mimicked by a 
static, spin-independent potential 



4.10 Spin splittings 

As for mesons, there are different possible approaches, and some compromises. 

Thanks to the charismatic role of Isgur and Karl, the chromomagnetic interaction is 
the best known mechanism to explain the spin-splittings of ground state baryons. It 
even works in sector where it should be prohibited. Take p — n: it gets a colour factor 
—16/3, a spin factor 4 (cr^c^ = +1 for triplet and —3 for singlet), and |</>(0)| 2 = (5^ 3 \r)) 
and results in about 0.6 GeV. For A — N, the colour factor is —8/3, the spin factor 6, and 
(5( 3 )(x)) is probably reduced, as the wave function is more extended. Hence a value 
A — iV of about 0.3 GeV is roughly compatible. 

The systematics of hyperfine splittings in a model based on 
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i<3 
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was remarkable. See, e.g., [114 139||180|181J . In particular, the chromomagnetic interac- 
tion, with its very specific (mimj)~ l dependence explains why the A is ligther than S, 
a long standing puzzle. See the paper by the Orsay group (Le Yaouanc et al.). This ap- 
proach has been regularly updated | 150|166| 182 - 185[ See, also, the interesting historical 
remarks by Lipkin [186], in particular on the pioneer contribution by A. Zakharov. 

The chromomagnetic model of hyperfine splittings has been used to make predic- 
tions for heavy baryons. It was for instance predicted that the £q — Aq splitting will 
remain large as the mass of the heavy quark increases, whil e Sq — Sq will decreases. 
See, e.g., 1 150 [, and for a recent upgrade in the beauty sector, 1 166 1. 

For the fine structure, the situation is less convincing. We have seen that the fine 
structure of mesons requires a sizeable spin-orbit interaction, which explains the order- 
ing 3 P < 3 Px < 3 P 2 . This spin-orbit potential is explained as to be resulting from a 
partial cancellation between the positive spin-orbit of the vector Coulomb potential and 
the negative spin-orbit of the scalar confinement. It has been argued [ 139[ 187fl that the 
cancellation might be more effective for baryons, but still the complete vanishing of the 
spin-orbit force is not understood. 
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There is definitely a problem with the spin-orbit splittings in baryons, which are 
usually very tiny, at least smaller than expected by extrapolation of meson spectroscopy, 
but sometimes surprisingly large. For instance the iV(1520) — iV(1535) difference (with 
J p = (1/2)- and J p = (3/2)-, respectively) is small, but (A(1520) - A(1405) is large. So, 
a strategy initiated by Isgur and Karl consists of setting to zero the spin-orbit forces in 
baryon, and to call for contributions beyond the simple quark model for the exceptions. 

Another approach, also present in the phenomenology of mesons, relies on the role 
of hadron loops. See, for instance, [188 189J . For sure, the A -h- nN coupling contributes 
to the A — N spacing. But the best advertised splitting is between A(1405) and A(1520). 
Both are negative-parity excitations of the ground-state A, with isospin 1 = 0. The 
A(1405) has J p = (1/2)-, and A(1520) has J p = (3/2)-. Hence, in the quark model, 
they correspond to the same spin structure as the ground state A, a light diquark with 
spin s q = giving a total quark spin s = 1/2. The motion of the s quark with respect to 
the (ud) diquark has i y = 1. Thus there are two possibilities, J = 1/2 and J = 3/2. The 
vicinity the KN threshold is presumably responsible for this unusually large separation 
A(1405) and A(1520). The pioneer here is Dalitz (190) . 

A third possibility consists of modifying the part of the quark dynamics dealing with 
the spin effects. For instance an instanton-induced interaction, inspired but the work of 
t'Hooft, was adapted by the Bonn group [138]. This was rather successful. Even more 
advertised is the class of models with a kind of pion-exchange between quarks inside 
the same hadron. This is named "Goldstone-boson exchange". 

Interesting models have been built were both chromomagnetism and Goldstone- 
boson exchange contribute to the spin splittings]^] 



4.11 Convexity properties 

There are not many exact results on the three-body problem. Everything becomes rather 
involved when going from two to three constituents. For instance, if you consider three 
distinguishable particles in a local potential, it is not obvious that the lowest energy for 
each spin J is an increasing function of Jp] 

Among the simple results, one may mention: the energies decrease if any mass is 
increased; . . . 

One point that was rather debated is the convexity as a function of the inverse 
masses. We have seen in the two-body problem that (the symbol (..) is meant for the 
ground state energy) that (m, m) + (M, M) < 2(m, M). Is is thus tempting to generalise 
as 

(m, m, rri) + (M, M, m!) < 2 (to, M, m') , (60) 



for any value of the "spectator" quark of mass vn! . This was conjectured [150 191 [. 



14 I remember a very interesting discussion with Glozman on this issue. Once a pion is emitted, the 
quark spin is flipped and thus the intrinsic spin-spin interaction at the quark level is modified. So it is 
perhaps not so easy to superpose within the same model Goldstone exchange and gluon exchange. 

15 I thank Andre Martin and Victor Mandelzweig for several discussions on this point. 
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After some investigation, it turns out that the rule remains true for reasonably-shaped 
potentials, but that there are counterexamples involving very sharp confinement and 
large mass ratios ||192 . 193[ . 

The lesson of ( [60) is that heavy quarks like to cluster together. 



4.12 Link between meson and baryon spectroscopy 
4.12.1 The baryon potential 

In the late 70s, there were clearly experts on baryon spectroscopy, obviously not re- 
stricted to heavy quarks, and on the other side, experts on quarkonia. But the question 
was already addressed of the link between the quark-antiquark interaction in mesons 
and the three-quark potential in baryons. 

One-gluon-exchange, or more generally, any colour-octet exchange gives a factor 1/2 
for the interaction of a pair of quarks in baryons as compared to the quark-antiquark 
interaction in mesons. This is the analogue of the (—1) factor for one-photon exchange 
in e + e + as compared to e + e~, or the Fermi-Yang rule G for the exchange of mesons, 
between iViV and NN , where G is the G-parity of the meson or set of mesons p*| 

It is thus reasonable to try, as an educated guess, to apply the same "1/2" rule to the 
whole potential. After all, when a quark in a baryon, say quark #1, is well separated 
from the two others, #2 and #3, it feels a colour 3 at separation r ~ r 12 ~ r 2 3, exactly 
as a quark feeling an antiquark at separation r, and the "1/2" rule reproduce perfectly 
this equivalence. 



This leads Stanley and Robson, Greenberg, Cohen, Lipkin, etc. |194-196[, to discuss 
the merits of this rule, and to propose a combined phenomenology of both meson and 
baryon sectors. I remember once Martin making an adventurous extrapolation of his 
(cc) and (bb) potential down to (ss) |64[. In the good old days of Inquisition, you would 



be burned at the stake for much less, because (ss) is anything but non-relativistic. Nev- 
ertheless, the agreement was good. I tried the "1/2" rule, and got the mass of Vt~(sss) 
within just a few MeV | |197[ . There are certainly many regularities in the hadron spec- 
trum that are mimicked by potential models, but the reason of not fully understood. 

The main issue is whether the interaction in baryons is pairwise. In nuclear physics, 
the interaction among nucleons is not pairwise: current potentials fitting the nucleon- 
nucleon scattering data and the properties of the deuteron underestimates the binding 
of few-nucleon systems. Three- and even four-nucleon forces have to be introduced. In 
atomic physics, there are empirical potentials between for instance two He atoms (Asiz, 
etc.). They can be reproduced by a (complicated but feasible) Born-Oppenheimer calcu- 
lation, where the energy of four electrons around two fixed a particles is properly min- 
imised. Then the interaction energy of three He atoms can be estimated, by minimising 
the energy of six electrons around three fixed a's. The interaction contains three-body 
terms, which are not very important in practice, but are there. 



16 The Fermi- Yang rule relates the potentials with the same isospin, for instance, pp and pn. 
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Similarly, any physical interpretation of the linear confinement in mesons leads to a 
generalisation for baryon which is of the three-body type. See, for instance, | |149 198 



|202[ , but many important papers have been probably omitted. In early days, some ap- 
proaches were a little empirical. For instance, in Jl49|, the same variant of the bag model 
used for quarkonia and hybrid mesons was applied to a bag containing three heavy 
quarks surrounded by their gluon field. At very large separations, the bag that min- 
imises the gluon energy is Y-shaped, with three tubes joining at the centre of the baryon. 
Today, the potential of baryons is derived in lattice QCD and even in AdS/ QCD, as will 
be explained by the other lecturers. 

The main result is that if you consider that V(r) ~ br in ( TT0| ) corresponds to a gluon 
tube of constant section whose linear behaviour oc r minimizes the energy, the analogue 
for a baryon is a F-shape set of three tubes whose cumulated length minimises the 
energy, i.e., 

V(r u r 2 , r 3 ) = b min(ru + r 2J + r 3J ) , (61) 



where the junction J is adjusted to minimise the sum. The solution is well known. For 
a flat triangle, with, say, Z213 > 120°, then J coincides with the quark #1, other wise, J 



sees each side at 120°, as shown in Fig. 16 




.4 



Figure 16: The string linking the quark to an antiquark in a meson (left) or the three 
quarks of a baryons has a minimal length. 



As a French lecturing in Italy, I cannot resist mentioning here the work of Fermat and 
Torricelli, and of Napoleon. The two first names will not be too much a surprise, and, 
indeed, Torricelli had an interesting correspondence with several scientists, including 
Fermat. Napoleon Bonaparte invaded Italy, among other countries, and his troops be- 
haved wildly, in particular in Naples, and the following sentence is reported "Non tutti 
i francesi sono ladri, ma buona parte si". Napoleon was also a fine mathematician, and 
his theorem on triangles is relevant for the confinement of baryons. 

The problem of the minimal cumulated distance to three given points was solved by 



Fermat and Torricelli. A possible construction is shown in Fig. 17 



One may also remark that the circumcircles of the external triangle such as BCA' 
intersect in J. The construction is also that of Napoleon's theorem, which states that the 



centres of these external triangles form an equilateral triangle. See Fig. 18 This is an 



interesting example of symmetry restoration starting from an asymmetric triangle. 
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Figure 17: If the quarks are located at A, B and C, and if A 1 is the third vertex of an 
external equilateral triangle built on the side BC, and similarly for B' and C , then the 
straight line AN , BB' and CC intersect at the Fermat-Torricelli junction. Moreover, the 
minimal sum of distances, Y = J A + JB + JC is equal to Y = AA' = BB' = CC 




A 



Figure 18: The Fermat-Torricelli junction is at the intersection of the circumcircles of the 
external equilateral triangles which are at the basis of Napoleon's theorem. 
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4.12.2 Inequalities between meson and baryon masses 

Presumably, the short range contributions are pairwise and fulfil the 1/2 rule, and the 
long range linear confinement is given by the Fermat-Torricelli minimal length Y (times 
the string tension). It is easy to check that the latter is bound by (p is the perimeter) 



V -<Y<^ 
2 " ~ 



(62) 



(the lower bound is saturated for a flat triangle, the upper one for an equilateral trian- 
gle), so that for the whole baryon potential, 



V > 2 K r i2) + v(r 23 ) + v(r 3 i)] . 



(63) 



In the 60s, physicists working on few-nucleon systems, or on quantum systems un- 
der gravitational interaction, have derived inequalities between three-body and two- 
body energies. In our case, 



V > -\v(r l2 ) + v(r 23 ) + v(r 31 )) 



and hence for the Hamiltonians 



2 m 



+ 



+ V > 



2 m 2 m V 12J 



+ 



(64) 



(65) 



Since the minimum of the sum is larger than the sum of minima, the above inequality 
immediately implies E 3 > (3/2) E 2/ or after adding the constituent masses |203[, 



2M(qqq) > 3M(qq) . 



(66) 



In other words, the matter is heavier per quark in the baryon form than in the meson 
form. The inequality ( |66] > has been generalised and improved in several ways. See, 
e.g., [204-206] for some details and references to the pioneer works, in particular by 



Hall, Post, Fisher, Ruelle, Levy-Leblond, etc. 

In the form (qqq) + (qqq) > 3 (qq), the inequality ([66]) implies that, at least within a 
simple quark model, baryon-antibaryon "annihilation" (more precisely rearrangement) 
into three mesons is energetically allowed, even at rest. 

On the other hand, if the mass ratio becomes large enough, the inequality is inverted, 
and becomes 

(QQQ) + (qqq)<3(Qq), (67) 

meaning that a triply-flavoured antibaryon would not annihilate into ordinary matter. 
The critical value of M/m is model dependent. There is in d67k another manifestation of 



the tendency of heavy constituents to cluster together, as seen already in p7[ ). 

"Why do it the easy way when you can do it the hard way?" was the favourite 
watchword of a popular TV series in the late 60s. An alternative (and not so rigorous) 
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way of understanding (qqq) + (qqq) > 3 (qq) is as follows. The reasoning will be useful 
in the section about multiquarks. 

Consider the case where the interaction is pairwise and the 1/2 rule is exact. Oth- 
erwise replace some equalities by inequalities. Both (qqq) + (qqq) and (qq) + (qq) + (qq) 
systems can be artificially considered as governed by a 6-body Hamiltonian (for more 
mathematical rigour, one can add a very weak overall harmonic confinement that makes 
all levels discrete and normalisable) 

H e [{9d}\ = E ^ + E Sid , E 9a = 3 • (68) 

i i<j i<j 

Let EQ[{gij}} denote the ground state. From the variational principle, it is maximal for 
the fully symmetric case where = 1/5 Vz, j. The baryon-antibaryon case corresponds 
to a set of coefficient ("*" means "times") {gij} = {6 * (1/2), 9 * 0}, and the three-meson 
configuration to {gj} = {3 * 1, 12 * 0}. The latter has the same mean value 1/5 as the 
former, but a much larger variance, about 0.17 vs. about 0.06, and this makes its energy 
lower. 



5 Multiquarks 

Le plus grand dereglement de Vesprit, 
c'est de croire les choses parce qu'on veut qu'elles soient, 
et non parce qu'on a vu qu'elles sont en effet. 

Bossuet 

5.1 Introduction 

We are now entering a delicate chapter. An atomic physicist could not feel satisfied with 
only the hydrogen atom and its variants, and the H and H 2 + ions, or a nuclear physicist 
with the deuteron, tritium and 3 He: they have to understand why some small systems 
are less easily bound than others, and to explore more complicated structures. Similarly, 
one has to know whether there exist hadrons beyond ordinary mesons and baryons. 
Besides multiquarks, several systems have been considered, either exotics, i.e., with 
quantum numbers that cannot be matched by (qq) nor (qqq), or crypto-exotics, with the 
same quantum numbers but different internal structure. The latter are more difficult to 
identify: one should first detect the exhaustive spectrum of ordinary hadrons, and look 
at super-numerous peaks, but the signature can be hidden by mixing of ordinary and 
extra-ordinary hadrons. 



17 The biggest disorder of the spirit, it is to believe things because we want that they are, and not because 
we saw that they are indeed. 
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5.1.1 Glueballs 



They used to be fashionable, and motivated several experiments. In a glueball, a few 
constituent gluons experience a confining interaction. There is a considerable phe- 
nomenology of glueballs, and many predictions based on lattice QCD, with noticeable 
fluctuations along the years |207[. 



5.1.2 Hybrid mesons 

At first sight, this is just another piece of the folklore, combining mesons and glueballs. 
However, it acquired some ground. 

Consider H 2 + in the Born-Oppenheimer limit, which is a very good approximation. 
At given pp separation r, one estimates the minimum of the electronic energy, which, 
when supplemented by the direct Coulomb repulsion, built the effective pp potential. 
In the Schrodinger equation, this potential gives the ground-state and the first orbital 
and vibrational excitations. If one takes the second electronic energy, one gets another 
effective potential generating a second set of levels. 

Now, replace the protons by a quark and an antiquark, and the electron by the 
gluon field, you first get the ordinary quarkonium and next the hybrids. Pioneer papers 
are | |208[|209[ . The quark-antiquark potential within an excited gluon field was explic- 
itly estimated in a variant of the bag model [ 210J , where gluons are treated classically 



and confined in volume whose shape and size is optimised: a QQ potential is derived, 
very similar to the Coulomb-plus-linear one, and a second potential is obtained, with 
a repulsive Coulomb interaction and a confinement smoother than for quarkonia. See 
Fig. [191 




Figure 19: Left: ordinary QQ potential, right: effective QQ potential for hybrids. 
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5.1.3 Hadron molecules 



It is regularly rediscovered that the Yukawa mechanism leading to the nuclear forces is 
by no means restricted to the nucleon-nucleon interaction. Any particle containing light 
quarks or antiquarks can emit or absorb light mesons and thus enter a nuclear-type of 
interaction. 

The hidden-charm sector is no exception. Already in 1975, Iwazaki suggested that 



ip' could be (ccqq) |211[ , where q denotes u or d. In 1976, Voloshin and Okun |212[ 



suggested the existence of molecules made of a charmed meson and an anticharmed 



meson. De Rujula, Georgi and Glashow |79[ were intrigued by one of the recurrence 
of J/ip, with the same quantum numbers J PC = 1 , whose branching ratios into DD, 
D*D + c.c. and D*D* were anamalous as compared to naive spin counting. They sug- 
gested a D*D* structure for the ^(4.03). In fact, the anomalies in the branching ratios 
were explained by the node structure of the state: roughly speaking, a two-body decay 
calls for a range of relative momentum, and the decay is either suppressed or enhanced 
if this momentum is near a zero or a maximum of the probability in momentum space. 

More recently, Tornqvist [ 213 214 [, and also, Ericson and Karl [215 |, Manohar and 



Wise 1 216 1, etc., stressed that the one-pion exchange between two charmed mesons or 
between a charmed and an anticharmed meson, can be attractive. This potential is 
somewhat weaker than the potential between a proton and a neutron, but of the same 
order of magnitude. It is experienced by heavy particles, and thus could lead to bind- 
ing. Remember that for an Hamiltonian H = p 2 /m + g V(r), where V(r) is attractive or 
contains an attractive part, the existence of a discrete spectrum depends on m g. 

The discovery of the X(3872) just at the D*D threshold was considered as the success 
of the prediction by Tornqvist and others. Some further measurements have confirmed 
the existence of this state. However, some decay properties, such as X — > ip(2S) + 7 
larger than X — > ip(lS) + 7, suggest the structure of a radially excited P-state. Pre- 
sumably both charmonium and molecular components exist in the wave function. The 
question now is whether this is just a charmonium state modified by the coupling to 
its nearest threshold, or whether there are really two separate states in the spectrum, 
one mainly 2 3 P 2 state of charmonium, another mainly a D*D + c.c. molecule, both with 
J PC = 

I cannot review the literature on molecules, especially the papers based on effective 
theories. A superficial reading gives the feeling that for any pair of hadrons, one can 
write up an effective Lagrangian that will lead to a molecule. The approach based on 
conventional nuclear forces is more predictive. Remember the Yukawa interaction in 
a S-wave contains a spin-isospin factor cr 1 .cr 2 Ti.r 2 , and thus is either attractive or re- 
pulsive and varies by a factor 1:9 in magnitude. After the discovery of the X(3872), the 
literature exploded. See, e.g., the reviews by Swanson | |217[ and by Nielsen et al. |218[. 



Recent contributions include [219 220 [ 
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5.2 Baryonium 



In the 70s, bumps were tentatively seen in antiproton cross sections, and peaks in inclu- 
sive 7-ray spectra of pp — > 7X. The generic name "baryonium" was invented for new 
mesons preferentially coupled to the baryon-antibaryon channels. For a review of the 



early results, see |221 [. 

A molecular picture of baryonium was proposed, e.g., by Shapiro and his collabora- 



tors 1 222 [, Dover et al. 1 223 - 225 [, etc. The interaction is given by meson-exchange (Fermi 



and Yang |226[ have given a simple rule to transform meson exchange into a particle- 



particle system to the corresponding particle-antiparticle system) at large distances, and 
exhibit a noticeable attraction in some partial wave. This suggests the possibility of NN 
bound states, in the same way that the exchange of mesons leads to a bound state of 
two nucleons. The role of annihilation remains however uncertain in this approach. 

A four-quark picture of baryonium was also elaborated, and in fact was even antici- 
pated by arguments based on the so-called "duality", whose explanation is beyond the 
scope of this review. See, e.g., Veneziano and Rossi | |221[ and Jaffe [ 227[ 228[. Unfortu- 
nately, if the authors of the above references were relatively careful in their predictions, 
this was not always the case with the followers, who produced dozens of baryonia and 
gener alis ations . 

Chan H.M., Hogaasen and their collaborators speculated about the role of the inter- 
nal colour degrees of freedom [229]J^] For (q 2 q 2 ), they named "true" baryonium a state 
where the quarks are in colour 3 and the anti quarks in colour 3. Such a state is thus 
easily coupled to the baryon-antibaryon channels but can be prevented of annihilating 
too fast if the quarks are separated from the antiquarks by some angular-momentum 
barrierj^] States with colour structure 66 were also envisaged and named "mock" bary- 
onia. It was suggested they could be rather narrow. Of course, what was not clear in 
this approach is why four-quark states should cluster in the form of a diquark and an 



antidiquark with such colour structure. See, e.g., [230]. 

No baryonium survived better experimental scrutiny, in particular with the intense 
and cooled antiproton beam of the LEAR facility at CERN. More recently, hadrons con- 
taining heavy quarks have been observed decaying into baryon + antibaryon + • • • , 
and enhancements are seen in the spectrum of the invariant mass pp, probably due to a 
strong final-state interaction. So, the saga of baryonium is perhaps not fully over. 

Note also that the mechanisms invented to explain the tentative baryonium states 
were applied to other configurations. For instance, Strottmann [ |231[ , Sorba and Ho- 
gaasen [232], and De Swart et al. (233J studied baryon states with four quarks and one 
antiquark (the name "pentaquark" was used later). Other studies dealt with (qq — qq) 
]234| or with dibaryons |235] . 



18 I thank Chan H.M., H. Hogaasen, P. Sorba and B. Nicolescu for numerous discussions 

19 Note that the structure (qq — qq) for angular momentum I > is not demonstrated, it is assumed 
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5.3 Chromomagnetic binding 
5.3.1 The chromomagnetic interaction 

As mentioned earlier, an appealing aspect of the quark model for mesons and baryons 



is the systematics of hyperfine effects such as J/ip — r\ c or A — N. Lipkin [ 186 236 [ 



likes to stress that Sakharov, though isolated, had interesting remarks on the subject. A 



popular reference is the paper of DGG |114[, where hyperfine forces among quarks are 



described in analogy with the Breit-Fermi interaction in atomic physics and attributed 
to one-gluon-exchange. See also ]115[ . 



This mechanism was used in the bag model \237\, and in several variants of potential 
models, with a component 

V ss = over Xi.Xj v ss (rij) , (69) 



rrii m, 
i<j ■> 

where v ss is either a contact term, to be treated at first order, or a regularised contact 
term, to be inserted in the wave equation. We already stressed the importance of the 
strength being proportional to (m; rrij)~ l , for instance for the E — A difference, but now 
we will first consider the case of equal masses. For angular momentum t > 0, there are 
additional terms, in particular a tensor component. 

5.3.2 The dibaryon H 



At the end of 1976, Jaffe |238| studied the (ssuudd) system in the SU(3)f limit, assum- 
ing that v ss = (v ss (rij)) has the same value for all the pairs, thus concentrating on the 
properties of the operator O = J2 (T i (T j K-^j i n this dibaryon system, as compared to 
the baryons. He discovered a very intriguing coherence, namely that for some configu- 
rations (O) can be attractive and larger (in absolute value) than its cumulated value for 
two baryons constituting the threshold! This is really remarkable. For the positronium 
molecule and for two positronium atoms constituting its threshold, both described by a 
potential Yl 9ij r i 3 l r we have the same cumulated strength Y gij = —2. 

Using for v ss the same value as in the light baryons, Jaffe predicted a "di-lambda", 
named H, bound by about 150 MeV below the baryon-baryon threshold. In the SU(3)f 
limit considered by Jaffe, the N, A, E or H receives from ( |69] > a downward shift of 
150 MeV (this is half of the A — N splitting), and thus the AA = NE thresholds ben- 
efit from about 300 MeV. Thus the H receives an additional 150 MeV as compared to this 
threshold. This was a dramatic change in the physics of multiquarks. Earlier specula- 
tions, such as baryonium, proposed metasable states, lying above their lowest thresh- 
old, but perhaps not too broad. The H, as initially proposed, was stable against any 
dissociation, and thus decaying weakly. 

As the MIT physicists are very convincing and powerful, the H was searched for in 
about 20 experiments]^] 

20 For comparison, the pentaquark analogue discussed in the next section was searched in a single 
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Unfortunately, it was later realized that the first calculation by Jaffe needed some 
improvements, namely | |239| - |242| : i) treat the spin-independent interaction (instead of 
assuming that (q 6 ) and (g 3 + g 3 ) starts from the same line before switching on the hy- 
perfine effects, ii) consider SU(3) F breaking that spoils the coherence of the configura- 
tions AA, NE, etc., , and ii) use for v ss values consistently computed from 6-body wave 
functions. Each effect dramatically reduces the binding, and the H becomes eventually 
unbound, though there is some model dependence. Recent lattice simulations [243] 244fl 
leaves the door open. 



5.3.3 The heavy pentaquark P 



In 1987, Gignoux et al. [245[ and Lipkin |246[ realised independently that the scenario 
for the H is repeated in a baryon configuration involving a heavy antiquark and four 
light quarks FT 

Let us consider a configuration (Qqqqq) in the limit mg for the antiquark and exact 
SU(3) p for the quarks. If the quarks have spin and form a triplet of flavour. This means 
(uuds), (udds) or (udss). Then the chromomagnetic operator O = (Ti-compo-j has 
an expectation value which is negative and twice larger than for A. This means that if 
you compare the pentaquark]^] P(Quudqs) to its thresholds (Qq) + (qqq), and adopt the 
same short-range correlation factor v ss , you find the pentaquark bound by aboundkarl 



150 MeV. This pentaquark was searched for in an experement at Fermilab 1 247 248 [, but 
the results were not conclusive. 

This pentaquark was further analysed, and its was found, as for H, that the cor- 
rections tend to spoil the binding. They are: breaking of SU(3)i?, finite mass for the 
heavy antiquark, self -consistent estimate of (v ss (rij)), or even better, inclusion of v ss in 



an accurate five-body calculation, etc. See, for instance, ]242 249 250 [ 



5.3.4 Other configurations 

The chromagnetic mechanism remains fascinating, even so it does not bind as easily 



as first thought. It contributes to the structure of the X, Y , Z states [251 [, and to the 
short-range hadron-hadron interaction. Leandri and Silvestre-Brac have scanned the 



configurations where interesting coherences occur [252 253J . Then they investigated 
whether the stability suggested by the chromomagnetic operator survived a more de- 
tailed calculation, with mostly negative results. 



experiment, and is vaguely in the agenda of the COMPASS experiment at CERN. "Selon que vous serez 
puissant ou miserable, les jugements de com vous rendront blanc ou noir. " (According as you're feeble, 
or have might, High courts condemn you to be black or white) Les animaux malades de la peste, Jean de 
la Fontaine 

21 Lipkin, in fact, had some precursor idea in this direction at Baryon 80 in Toronto. Anyhow, I remem- 
ber Lipkin arriving at Les Houches for the Workshop on the "Elementary Structure of Matter", handing 
me a freshly issued preprint, myself welcoming him with a copy of our latest manuscript, and our aston- 
ishment when we realized that both papers were dealing with the same exotic baryon! 

22 The word was adopted at this occasion 
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5.4 Binding by central forces 

In a novel by Alphonse Daudet, one finds this sentence: " . . . il ferma la porte a double 
tour. Malheureusement, il avait oublie la fenetre ..." (He double locked the door, but 
forgot the windows). Similarly, if the chromomagnetism cannot bind multiquarks, why 
not try chromoelectricity? 

Believe it or not, it works, and there is an amusing side effect, or I should perhaps say, 
a miracle: while the literature on multiquarks is usually a jungle, of controversed, here a 
consensus has been reached that some configurations should support stable tetraquarks, 
i.e., states below the threshold for their dissociation into two mesons. 

I will first summarise the results obtained from an empirical pair- wise ansatz for the 
multiquark interaction, and in the next subsection, explain the changes brought by an 
improved modelling of the dynamics. 



5.4.1 The additive model 



As explained in Sec. 4.12 the potential for mesons can be extrapolated to baryons by 
the "1/2" rule, i.e., V = [v(r 12 ) + v (r^) + v (r 3 i)]/2, where v is the quarkonium potential. 
It can be read as the exchange of a colour octet in the t-channel. Similarly, in nuclear 
physics, there are isoscalar exchanges, leading to an isospin-independent interaction, 
and isovector exchanges, which get a ti.t 2 factor, such a pion-exchange. The analogy 
is of course delicate, as local and global symmetries are not on the same footing. Nev- 
ertheless, the multiquark phenomenology in simple constituent models has been very 
often based on the assumption that 

V = -—Y;X i .\ j v(rv), ( 7 °) 

which reduces to the 1/2 rule in the case of baryons. Here is the version for quarks. For 
an antiquark, \ — > —A* 

This additive model leads to coupled equations, using either a (qq) — (qq) basis with 
colour 1 or 8 for the pairs, or a (qq) — (qq) basis with internal colour 3 — 3 or 6 — 6. 

Solving the four-body problem is a considerable task, and requires a lot of care. It 
can be done by variational methods based, for instance on a superposition of Gaussian 
wave functions, or by hyperspherical expansion. 

The hyperspherical expansion consists of considering the n Jacobi variables pjof a 
(n + l)-body problem, as a unique 3 n-dimensional vector, and to use spherical coordi- 
nates (q, 0) to describe it. In the particular case of the harmonic oscillator with equal 
masses, the potential is isotropic, V oc Yl r ij x Q 2 > an d a single radial equation allows 
one to recover the standard results for this interaction. In general, one deals with an 
anisotropic potential in the 3 n-dimensional space, and thus an infinite set of coupled 
equations (the same situation would occur for a single paticle in a potential that is not 



23j f 

necessary rescaled to be associated with the same reduced mass 
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rotation-invariant). One obtains a rather satisfactory convergence as a function of the 
number of equations actually included. 



5.4.2 Results for equal masses 

If this colour-additive potential, without spin-dependent terms, is applied to (q, q q, q), 
no bound state is found below the (qq) + (qq) threshold. Even though such a model is 
rather primitive, it already indicates that there is no proliferation of narrow multiquarks. 

One may wonder why the positronium molecule is stable whereas the (q, q, q, q) is 
not in such a simple quark model. The generic Hamiltonian for these systems is 

^— ' 2m 8m / — ' z — ' 

l Kj Kj 

where v(r) — — 1/r for the atomic problem, and v(r) = \ r + ■ ■ ■ in the quark problem. 
For the thresholds (atom-atom or meson-meson), the molecule, and the "true" T = (33) 
and "mock" M = (66) colour wave-functions of the tetraquark system, the coefficients 
are given in Table |4j 



Table 4: Coefficients of the potential energy for the threshold, the positronium 
molecule and the T and M type of tetraquark, average and variance 



Sta*>^^ 


12 


34 


13 


24 


14 


23 


9 


A<7 


Threshold 








1 


1 








1/3 


0.22 


Ps 2 


-1 


-1 


1 


1 


-1 


-1 


1/3 


0.89 


T 


1/2 


1/2 


1/4 


1/4 


1/4 


1/4 


1/3 


0.01 


M 


-1/4 


-1/4 


5/8 


5/8 


5/8 


5/8 


1/3 


0.17 



We repeat here in some more detail a reasoning sketched at the end of Sec. 4.12.2 
dealing with the comparison of mesons and baryons. Symmetry breaking is known to 
lower the energy of the ground state. For instance, the linear oscillator h(0) = p 2 + x 2 
has a ground state at e(0) = 1 which becomes e(A) = 1 — A 2 /4 for h(X) = h(0) + Ax. 
More generally, for H(X) = H + A Hi, where H is even and Hi odd, one gets the same 
lowering e(A) < e(0) of the ground state, and this is true for symmetries other than 
parity. This is easily seen by using the variational principle, with the symmetric ground 
state of H as a trial wave function for H(X). 



Thus, if one considers the family of Hamiltonians ( 71 ), the ground state energy is 
maximal if all couplings are equal, i.e., <7jj —g and Ag = 0. 

If Hi = H + X Hi with Hq symmetric and Hi a given asymmetric potential energy, 
then the ground state e(A) is maximum at A = and concave. Thus the larger A for 
A > 0, the lower the energy. And the larger | A | for A < 0, the lower the energy. 



MULTIQUARKS 



57 



This is rigorous, so far. But the energy, as a function of A is nearly symmetric, and for 
H 4 = H +X Hi, the energy e(A) decreases approximately as e — A 2 e . Furthermore, in the 
case of ( f7T) , H corresponds to the symmetric case g^ = g, and one can probe different 
Hi, i.e., different ways of breaking the symmetry among the pairs in the potential: the 
energy decreases approximately according to the variance Ag of the distribution of the 
coefficients {gij}. 

From the values listed in Table |4| one sees that Ag is larger for Ps 2 than for two non 
interacting atoms. This explains (again, not very rigorously) why PS2 is stable in atomic 
physics. On the other hand, for the T or M tetraquarks, Ag is smaller than for two 
non-interacting mesons. This suggests that in the simple additive model without spin- 
dependent terms, there is no stable tetraquark, and the mixing of colour configurations 
does not improve much the situation. In other word, the tetraquark is penalised as 
compared to the positronium molecule by the non-Abelian character of the theory. 



5.4.3 Binding for unequal masses 

The lesson of the above subsection is that the tetraquark potential is too symmetric 
with respect to permutations, as compared to the potential governing the two-meson 
threshold. So to get binding, another effect is required. We have mentioned chro- 
momagnetism, we shall see below the possibility of improving the modelling of spin- 
independent forces, to better conform with QCD. There is another effect, which is con- 
ceptually simpler: introduce symmetry breaking in the kinetic energy, i.e., mix different 
flavours in the multiquarks. 

As stated earlier, any symmetry breaking lowers the ground-state energy, but does 
not necessarily improve binding. Usually the threshold benefits more from symmetry 
breaking, thus the stability deteriorates. 

This is the case, for instance, in atomic physics, when one breaks permutation sym- 
metry for For simplicity, let us do it in the same way in the two charge 
sectors, i.e., keeping an eigenstate of charge conjugation. This corresponds to (M + , m + , 
M, m~), and one can use the scaling properties of the Coulomb problem to impose 
M^ 1 + vrC 1 = 2 ^r 1 without any loss of generality. This corresponds to maintaining 
the threshold energy constant. What is observed is that the positronium-line molecule 
is stable for M = m, but becomes unstable for M > 2.2m or m > 2.2 M. For a large 
value of the mass ratio, the ground-state consists of a protonium atom which cannot 
polarise a positronium atom to attach it. (Of course, some metastability of an hydrogen 
atom and an antihydrogen atom cannot be exclude, but this is a much more delicate 
issue.) 

The same mechanism renders (Q,Q,q,q) less bound than the equal-mass tetraquark 
in any simple, flavour-independent, central potential. This is of course rather unfor- 
tunate, since the best candidates today, to be discussed in Patricia's lectures following 
mine, are just in the hidden charm or hidden beauty sector! The mechanism here is 
presumably long-range forces, such as pion exchange, and spin-dependent forces. 

Another possibility if symmetry breaking with two masses consists to give up charge 
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conjugation and keep particle identity in the two sectors. In atomic physics, the config- 
uration reads (M + , M + , m~,m~). It is well known that the binding energy of the hydro- 
gen molecule (in units of the threshold energy to have a scale-independent statement) is 
much larger than that of the positronium molecule. Historically, the hydrogen molecule 
was studied in the Born-Oppenheimer-Heitler-London limit, starting from M — > oo, 
while the question of the stability of the positronium was addressed later by Whyleer, 
Ore, Hylleraas, etc. 

The Hamiltonian can be written as follows 

H(M + ,M + ,m-,m-) = H( l x + ,^, l x-^-) + ^(^-j^[pj+pi-pl-pl] , (72) 



In ( 72 1, with the choice of the inverse of /i averaging the inverse masses, the two Hamil- 
tonians have the same threshold, but the hydrogen molecule benefits from the term 
breaking charge-conjugation, and thus has a lower ground state! In other words, the 
stability of the hydrogen molecule is a consequence of the stability of the positronium 
molecule, via the variational principle. 

However, the Coulomb character of the potential hardly matters, what is crucial is 
property that the same potential governs both systems. Thus it is not a surprise that 
in a simple quark model with flavour independence, one finds that (Q,Q,q,q) becomes 



stable if the mass ratio M/m is large enough. This was found by Ader et al. 1 203 [ and in 
many subsequent papers, see, e.g., [ 220| 254 - 261J , and references there. 



The question is now how large M/m should be. It depends on the details of the 
model, it depends how good is the four-body calculation, and whether or not spin- 
dependent corrections are also added. For long, it was believed that with current pair- 
wise models, (6, b, q, q) was required, where q — u, d. Two b or not two b? The question 



was answered by Jang and Rosina 1 259 [, who got binding with Q = c, using an additive 



potential model fitting mesons and baryons. This was confirmed by Barnea et al. |261 [. 



5.5 Multiquarks in a Steiner-tree potential 
5.5.1 A new confinement for tetraquarks 

As already discussed, the pairwise interaction with colour factors is not dictated by 
QCD. It is a convenient tool for preliminary investigations. For sure, the ansatz looks 
reasonable for a iV-body colour-singlet, when a quark is well separated from the (N — 1) 
remaining constituents: the quark feels the same interaction as in a quarkonium. But 
there is no fundamental reason why the interaction should be pairwise. 

Years ago, Artru p98| , followed by many others pl^[l99f|202|[262| , pictured the 



mesons as a simple string linking the quark to the antiquark, and baryons as three 
strings, each linking a quark to an intersection named junction. In the strong coupling 
regime of QCD, this corresponds to a flux of gluons linking the quarks. The flux having 
a constant cross section, the potential grows linearly. 
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The analogue for tetraquarks of the F-shape interaction of baryons consists of two 
terms, 



1. A flip-flop interaction, Vf, which is the minimum of v (ri 3 ) + v (r 24 ) and v(ru) + 
v(r 23 ), 

2. A connected Steiner-tree, V s , which links the quarks to the antiquarks with a min- 
imal length. 



The combination of flip-flop and Steiner tree, with a linear behaviour, namely, 

V^min^K) , 



(73) 



thus involving a discrete two discrete minimisations and a minimisation over the con- 
tinuous variables describing the location of the junction of the Steiner tree. This inter- 
action was proposed by Carlson et al. |263[, Vijande et al. |264[, and Ay et al. [265 1. The 
model ( [73] > is schematically pictured in Fig. 20 



X 



Figure 20: Simple string model for the tetraquark: the potential is the minimum of two 
meson-meson links (flip-flop) and a connected contribution. 



5.5.2 The flip-flop interaction 

The flip-flop interaction was proposed in a celebrated paper by Lenz et al. (2661, with, 
however, a quadratic behaviour for v(r). In this paper is discussed the possibility of 
extending the quark model of ordinary mesons and baryons to the study of multiquark 
states and of the hadron-hadron interaction. 



5.5.3 The connected string 

It has interesting properties. In the planar case, it can be constructed by duplicating 



Napoleon's construction for baryons, see Fig. 21 



In space, the length of the minimal Steiner tree is the maximal distances between 
two circles. The first one includes all the points making an equilateral triangle with the 
quarks, the second one for the antiquarks. See Fig. 22 Now if you wish for the best algo- 
rithm to estimate the maximal distance between two circles in space, you have to read 
the reviews on computer-assisted cartoons. The software instructions in this branch 
of art are such that once a few basic elements are drawn, you are guided to deduce 
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Figure 21: In a plane, the Steiner-tree potential can be constructed as the Fermat- 
Torricelli minimal sum of distances for baryons. The length of interest is equal to dis- 
tance between the two auxiliary points u>i 2 and w 3i . More precisely, to the quarks v± 
and i> 2 are associated the Melznak set s = {wi 2 , h 2 } and to the antiquarks v 3 and v 4 are 
associated s = {w 34:1 1 34 }, so that the triangles (u>i 2 , i>i, v 2 ) and (1034, i>3, t^) are equilateral. 
The length of the minimal Steiner tree linking v±, v 2 to t> 3 , t> 4 , which involves the opti- 
misation of several continuous variables is equal to the maximal distance from s to s, 
which is just a discrete optimisation ,and thus very much faster. 



auxiliary pieces of the picture. So experts on multiquark physics can seek positions at 
Hollywood! 

A simplified version of the connected interaction V s can be found in the work by 
Barbour and Ponting [267,268]. A set of Jacobi coordinates is shown in Fig. 23 It reads 
(besides the centre of mass) 



x = r 2 - r 1 



y 



r 4 - r 3 



r<i + r 3 -r 2 - n 
x/2 



(74) 



so that the intrinsic part of the four-body Hamiltonian reads 

,2 „2 9 



m m 



PI 

m 



V . 



(75) 



In the case of a T-baryonium, with colour 3 coupling of the two quarks, and 3 for the 
antiquarks, Barbour and Ponting adopted a confining interaction 



V = a(x + y + z/y/2) . 



(76) 



Then H A exactly splits into three independent pieces, and the ground state is at e a 



2/3 



m 



-1/3 



[2 + 2 x / 3 ] above the threshold at 2 eo a 2 / 3 m 1 ^ 3 . Here eo is the negative of the first 



zero of the Airy function, see Sec. 3.3 



Now, if one gives up the square geometry of Barbour and Ponting, and minimises 
the string potential by varying the location of the junctions, one gets a lower energy, but 
still above the threshold, for equal masses. 
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Figure 22: The minimal length of the Steiner tree viSi + v 2 S\ + Sis 2 + s 2 v 3 + s 2 t>4 linking 
the quarks to the antiquarks is the maximal distance between their Melznak's 

circles. Any point of the circle C\ 2 makes an equilateral triangle with the quarks v 1 and 
v 2 . The same holds in the antiquark sector. 



'7 



'I 



Figure 23: Confining interaction for a tetraquark in the model by Barbour and Ponting : 

V = a (x + y + z) 
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5.5.4 Rigorous results 

We now discuss the four-body energy with the full linear potential V 4 of Eq. ( |73] >, starting 
with the exact result, and then, in the next subsection the numerical results. 

There are not many exact results for the four-body problem, and with the cumber- 
some potential ((73]), it looks difficult to derive a simple upper bound. In particular, V4 is 
smaller than the flip-flop Vj alone, which in turn is smaller than the potential a ri 3 + a r 24 
which governs the threshold. This means that the effective meson-meson interaction is 
attractive. In a world with on or two dimensions, this would imply binding. But in our 
valley of tears in three dimensions, this is not the case. 



Nevertheless, it is easily found from the construction of Fig. 22 that 



V s <a 



{x + y)— + zV2 



(77) 



and it has been shown [265 1 that this bound also holds for V4. Thus one gets a upper 
bound 



Ha < 



m 



m m 



+ a 



(78) 



which is exactly a sum of three independent terms. Is is easily found that in this model 
(QQqq) becomes stable for a mass ratio M/m > 6402. This is of course very crude, as the 



inequality in the potential is saturated only when the circles in Fig. 22 are in the same 
plane, and do not overlap. 



5.5.5 Numerical results 

The four-body problem with the interaction Va can be solved using standard techniques. 
Here, the difficult y comes from the numerical estimate of the matrix elements of the 
potential. It is found that binding is obtained even for equal masses |264|, and, of course, 
it becomes better for (QQqq). However, the minimisation in ( |73~| and in the flip-flop 
itself implies a continuous rotation of the internal colour degree of freedom. Hence 
colour cannot contribute to the antisymmetrisation. In other words, the stability if for 
distinguishable quarks and antiquarks. The connected term, or Steiner tree, contributes 
marginally, as compared to the flip-flop. 



5.5.6 String model of pentaquark 

The model can be adapted to the case of pentaquarks, i.e., four quarks and one antiquark 
1 269 1 . The potential is pictured in Fig. 24 Again, the interaction contains flip-flop terms, 
i.e., the lowest of the meson + baryon string potential, and some connected diagrams, 
the rule being to take the minimum. 

Again, the calculation corresponds to distinguishable quarks, even in the case where 
for simplicity the same mass is adopted. Then it is found that the pentaquark is stable. 
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Figure 24: Simple string model for the pentaquark configurations 



5.5.7 String model of hexaquarks 



1270]. 



In a recent study, six-quark configurations were studied, both (q 3 q 3 ) and (g 6 

The string potential for (q 3 q 3 ) is shown in Fig. 25 Missing is a term consisting of 



a meson and a separated tetraquark. For quarks and antiquarks with equal masses, 
but distinguishable, a weak binding is obtained. The lowest threshold is (qq) 3 . If one 
introduces two different masses and look at (Q 3 q 3 ), the binding first is improved when 
M/m increases from M/m = 1. But for very large M/m, the lowest threshold becomes 
(QQQ) + (qqq) and (Q 3 q 3 ) becomes unstable. 






Figure 25: Potential for (q 3 q 3 ): disconnected terms of baryon-antibaryon type of three- 
meson type (with suitable minimisation over the permutations) and some connected 
terms. 



In the case of (q e ) the potential is either of baryon-baryon type, with all possible per- 
mutations, or given by a connected string. See Fig. 26 For equal masses, a comfortable 
binding is obtained. For two different masses, i.e., (Q 3 q 3 ), the binding deteriorates. As 
the mass ratio M/m further increases, the 6-quark system becomes unbound and the 
ground state consists of two isolated baryons, (QQQ) + (qqq)- 



5.5.8 The next challenge: combining string and full antisymmetrisation 

The above string model, with mainly a flip-flop interaction and some connected contri- 
butions, has many appealing features. However, it is not suited for implementing the 
constraints of antisymmetrisation, and it tends to produce several stable multiquarks 
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Figure 26: Contributions to the linear part of the (g 6 ) potential, disconnected and 
conected terms. A minimisation over the permutations is implied. 

that are not physical. Consider for instance the tetraquark case. The potential is taken as 
the minimum (in units of the string tension) of r 13 + r 2 4 which corresponds to a singlet- 
singlet state of colour, or r 14 + r 2 3 which corresponds to another singlet-singlet, and the 
double-F string which corresponds to (3 — 3)i. For distinguishable quarks and anti- 
quarks, this is plausible, and the string potential is just a Born-Oppenheimer potential, 
once the coloured gluon field is minimised. Some change is required to implement this 
dynamics with full account of the symmetry constraints. 



5.6 The light pentaquark 



The claim by the LEPS collaboration |271[ that there exists a rather narrow baryon 9 + 
with mass about 1.54 GeV, charge Q = +1 and strangeness S = +1 stimulated much 
activity in the physics of hadron. 

First, other experimental groups looked back at their old data already on tapes, and 
also found some evidence for this state or for its SU(3) p partners (In the paper by Di- 
akonov et al. |272[ having motivated the experimental search, the 9 + supposedly be- 



longs to an antidecuplet 10). Later, the experimental evidence became much weaker. 
For a summary, see, e.g., the review by PDG []. 

Second, the original calculation within chiral dynamics was revisited 1 273 2741 and 
led to different conclusions. The pentaquark was also studied with QCD sum rules or 
Lattice QCD, with a variety of conclusions. There were also studies within potential 



models, with detailed five-body calculation. See, for instance, |275[. The light pen- 
taquark can not be reproduced in such calculations, unless one introduces ad-hoc extra 
assumptions. 
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6 Outlook 

Important problems never receive a definite solution 

Confucius 

About 50 years after the first speculations on quarks, what does remain? Let us just list 
of few items for the discussions after the lectures. 

SU(3) F symmetry is now becoming somewhat controversial. Sometimes, the SU(3) F 
symmetry is almost exact. This is the case for the rates of J/ip — > pp, AA, .... In other 
cases, SU(3)f looks badly broken. For instance, if one looks at K to 7r production in high- 
energy collisions, one finds a large "strangeness suppression". Perhaps the elementary 
couplings for uu or ss production are very similar, but the small differences in masses 
are amplified by a tunnelling factor which is exponential, like in the Gamow theory of 
a decay. 

From flavour symmetry to flavour independence In hadron spectroscopy, it is del- 
icate to attempt a blind extension of SU(3)f to SU(4) or SU(5) to include charm and 
beauty. The breaking is too large to be treated as a perturbation. A new understanding 
of the symmetry has been provided by the quark model. The basic interaction is flavour 
independent, or has a well-defined flavour dependence in the spin-spin and spin-orbit 
terms. Then the same basic interaction is used to build the various (qq') mesons and 
(qq'q") baryons. 

The bootstrap idea was conceptually superb, but operationally somewhat inefficient. 
Still, a kind of duality exists^} for instance, it is natural to consider the A baryon in 
some context the (qqq) partner of A^ = (qqq) with all quark spins aligned, but to describe 
the pion-nucleus scattering, there is nothing better than A as a n — N resonance. When 
the quark model was invented, it was considered as a relief, with at last a systematic 
and coherent scheme to describe the hadrons. Nowadays, for any hadron a of any mass 
m, one finds a pair of existing hadrons b and c and an effective Lagrangian coupling a, 
b and c, such that a appears as a kind of molecule. The overall picture is seemingly lost. 
Is that a big step backwards? 

From quark models to QCD This will be covered in some of the next lectures. It is 
fascinating how a somewhat empirical model, but based on the acute observation of 
facts, became a beautiful theory. 

A premium to simplicity The quark model, even in its simplest variant, works rather 
well. Presumably, the quark model incorporates a lot of subtle dynamical effects in an 
effective way, and enables one to make successful predictions. 

24 G.F. Chew, private communication, years ago, that is very thankfully remembered 
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Hadron dynamics beyond hadron decay An ambitious application of the quark model 
deals with the hadron-hadron interaction. This is the continuum part of the sections de- 
voted to loop corrections and to multiquark spectroscopy. 

Let us mention the calculation of the KN interaction in the exotic S — +1 sector 
of strangeness, and the extensive work done to understand the short-range part of the 
nucleon-nucleon interaction in terms of quarks. 
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